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ABSTRACT 

We investigate the contribution of star-forming galaxies to the ionizing background at z ~ 3, building 
on previous work based on narrowband (NB3640) imaging in the SSA22a held. We use new Keck/LRIS 
spectra of Lyman break galaxies (LBGs) and narrowband-selected Lya emitters (LAEs) to measure 
redshifts for 16 LBGs and 87 LAEs at z > 3.055, such that our NB3640 imaging probes the Lyman- 
continuum (LyC) region. When we include the existing set of spectroscopically-conhrmed LBGs, our 
total sample with z > 3.055 consists of 41 LBGs and 91 LAEs, of which nine LBGs and 20 LAEs are 
detected in our NB3640 image. With our combined imaging and spectroscopic data sets, we critically 
investigate the origin of NB3640 emission for detected LBGs and LAEs. We remove from our samples 
3 LBGs and 3 LAEs with spectroscopic evidence of contamination of their NB3640 flux by foreground 
galaxies, and statistically model the effects of additional, unidentified foreground contaminants. The 
resulting contamination and LyC-detection rates, respectively, are 62 ± 13% and 8 ± 3% for our LBG 
sample, and 47 ± 10% and 12 ± 2% for our LAE sample. The corresponding ratios of non-ionizing UV 
to LyC flux-density, corrected for intergalactic medium (IGM) attenuation, are 18.0t=5 4 4 8 for LBGs, 

and 3.7tH for LAEs - We 

use these ratios to estimate the total contribution of star-forming galaxies 
to the ionizing background and the hydrogen photoionization rate in the IGM, finding values larger 
than, but consistent with, those measured in the Lya forest. Finally, the measured UV to LyC flux- 
density ratios imply model-dependent LyC escape fractions of f^ s y c ~ 5 — 7% for our LBG sample 
and f/f L c - 10 - 30% for our fainter LAE sample. 

Subject headings: galaxies: high-redshift - intergalactic medium - cosmology: observations - diffuse 
radiation 



1. INTRODUCTION 

Identifying the sources of the radiation that re-ionized 
the intergalactic medium (IGM) at redshift z > 6 is a 
key challenge for observational cosmology. Massive stars 
in star-forming galaxies are considered the most likely 
source of the Lyman-continuum (LyC) photons necessary 
for reionization. QSOs, while also efficient producers of 
ionizing radiation, appear to be too rare at high redshift 
to produce suffici ent ionizing flux (jHopkins et al.l [20071 
ICowie et al.l 120091) . Unfortunately, direct observation of 
ionizing radiation from high-redshift star-forming galax- 
ies is not possible, as the Universe remains opaque to LyC 
photons above redshifts z ~ 4 due to the abundance of 
residual neutral gas in the IGM. Thus, empirical con- 
straints on the sources responsible for reionization come 
primarily from (i) the determination of the non-ionizing 
UV luminosity function of the highest-redshift galaxies 

1 Based, in part, on data obtained at the W.M. Keck Obser- 
vatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California, 
and NASA, and was made possible by the generous financial sup- 
port of the W.M. Keck Foundation. 

2 David and Lucile Packard Fellow 



currently observable, and (ii) direct measurements of ion- 
izing radiation from lower-redshift analogs of the z > 6 
galaxies responsible for reionization. 

Progress has recently been made in both of the 
above approaches. Bright, color-selected z ~ 6 - 7 
candidate galaxi es have recently been spectroscopically 
confirmed fe.g IPentericci et aLI l20lTl: IQno et al.l [20121: 
ISchenker et al.ll2012h . Searches for fainter, line-emitting 
galaxies using narrowband imaging techniques have un- 
covered significant samples of Lya e mitters (LAEs) at 
redshifts as high as z ~ 7 (e .g., IKrug et ail 120121 : 
IQta et alj|20iq lOuchi et al.ll2010D . At even higher red- 
shift, deep near-IR observations of the Hubble Ultra 
Deep Field using the WFC3 camera on the Hubble Space 
Telescope (HST), combined with optical H ST /ACS 
data, have re vealed populations of redshift z ~ 7 — 10 
galaxies (e.g.. IBouwens et al.l 120101 , IBunker et al.l 120101: 
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Among the results of such studies, the emerging 
z ~ 8 luminosity functions suggest that there are too 
few bright galaxies to reionize the IGM, implying that 
relatively low-mass galaxies are important contributors 
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to the ionizing flux budget. 

The effort to unambiguously identify low-redshift 
analogs of the galaxies responsible for reionization has 
been difficult. Observations to directly detect ionizing 
continuum escapin g from galaxies at < z < 2 have been 
unsuccessful (e.g.. jMalkan et al.ll20Cjl iCowie et al Il2009t 
iGrimes et a!]|2009HSiana et al.ll2010l: iBridee et afl koiO). 
Thus, searches have turned to the highest redshifts at 
which the opacity to LyC photons through the IGM, 
T LyC, is is 1) i- e -: z ~ 3. The expected rest-frame UV 
flux levels are low at such redshifts, requiring deep ob- 
servations to identify sources and to detect escaping LyC 
emission. The low flux levels also complicate the inter- 
pretation of presumed detections of LyC flux, as the sky 
surface density of faint foreground sources is large (~ 75 
arcmin -1 ), leading to the possibility of contamination 
from lower-redshift sources. With these caveats in mind, 
there have been several reported detections of escaping 
LyC flux from galaxies at z ~ 3 throug h both spec- 
trosc o py of Lyman break galaxies (LBGs) ([Steidel et al.1 
2001: iShaplev et all I2006D and narrowband imaging of 
LBGs and LAEs (llwata et all [20091: iNestor et all 120111: 
IVanzella et al.H201lD. 

In INestor et al.l ( 20111) . we searched for escaping LyC 
flux from z > 3.055 LBGs and photometrically-selected 
LAE candidates in the SSA22a field, whi ch contains a 
large over-density of galaxies at z = 3.09 (jSteidel et al.l 
1998). In that paper we reported the detection of six z > 
3.09 LBGs and 27 candidate z ~ 3.09 LAEs through a 
narrowband filter that is opaque to non-ionizing flux from 
sources at z > 3.055, and thus probes light blueward of 
the Lyman limit for sources at z ~ 3.09. We interpreted 
these detections as direct evidence of escaping LyC flux, 
which we in turn used to estimate the comoving ionizing 
emissivity, CLyCt a t z ~ 3. Our primary conclusions were 
that the contribution to tLyC from star- forming galax- 
ies at z ~ 3 exceeds, but is roughly consistent with, 
given our uncertainty, that expected from determina- 
tions of the photoionizatiqn rate in the Lya forest (e.g . , 
iBolton fc Haehnelti [20071 iFaucher-Giguere et al.l 120081 ). 
and that the bulk of this contribution comes from rel- 
atively faint sources (Mab ^ —20) such as those that 
compose our sample of LAEs. 

Our ability to constrain tLyCi however, was limited 
by our fairly small sample of LBG detections, which is 
very likely to cont ain contamination b y foreground in- 
terlopers (see, e.g.. IVanzella et aL[|2012fl . As our seeing- 
limited observations restricted our ability to account for 
contamination by foreground galaxies in individual sys- 
tems, we used a statistical approach for th e sample as a 
whole . Furthermore, many of the results in INestor et al.l 
(2011) were based on the sample of LAEs, the majority 
of which were not spectroscopically confirmed. Each of 
these limitations can be alleviated with follow up spec- 
troscopy. Therefore, in order to increase the size of our 
LBG sample, we obtained spectra of additional color se- 
lected LBG candidates lacking prior spectroscopic red- 
shift confirmation. We also obtained spectra of as many 
of our narrowband-selected LAE candidates as possible, 
in order to confirm their redshifts. These data have the 
additional benefit that, in most cases, we can directly 
search for evidence of foreground contamination in in- 
dividual objects. We also obtained particularly deep 



spectra of many of the sources having apparent LyC de- 
tections, which allow us to perform a detailed analysis 
(e.g., of the spatial distribution) of the detected Lya, 
non-ionizing UV, and LyC fluxes. 

The organization of this paper is as follows. In Sec- 
tion d we describe the observations and reduction of the 
data used in this study. We update our LBG and LAE 
samples based on the results of our new spectroscopy in 
Section [3] In Section [4] we discuss individual LBGs and 
LAEs with presumed LyC detections. In Section [5] we 
describe our techniques for statistically accounting for 
foreground contamination and attenuation of LyC flux 
by the IGM. We present updated estimates of CL y c and 
discuss LyC escape fractions in Section[6l and summarize 
our results in Section [7] Throughout the paper all magni- 
tudes are in the AB system, and we assume a cosmology 
with fl m = 0.3, fl A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Photometric Observations and Sample 

Our analyses make use of deep multiband imaging 
available in the SS A22a field. These da ta, which are de- 
scribed in detail in INestor et al.l ()201lD . include ground- 
based broad B-, V- and i?-band images and narrow- 
band images with effective wavelengths at A ~3640A and 
A -4980A (hereafter NB3640 and NB4980, respectively). 
The NB3640 and part of the NB4980 data were obtained 
using the Keck/LRIS imaging spectrograph, while the 
broadband imaging and the remainder of the NB4980 
data were obtained with the Subaru/Suprime-Cam. Ad- 
ditionally, archival H ST /ACS F814W imaging is avail- 
able for 70% of the field. 

The NB3640 filter has a central wavelength of 3635A 
and FWHM of lOOA. For sources at z ~ 3.09, NB3640 
samples the rest-frame spectral range A ~ 875 — 900A. 
It is opaque to wavelengths longward of the redshifted 
Lyman limit for sources above z ~ 3.055 and therefore 
provides a clean probe of escaping LyC emission for such 
galaxies. The effective wavelength of the broad i?-band 
filter (A ~ 6510A) corresponds to rest-frame A ~ 1600A 
at z ~ 3.09. The NB3640— R color is accordingly a 
measure of the non-ionizing to ionizing UV flux-density 
ratio, Fuv/FLyC, for galaxies with z > 3.055. The 
ACS/WFC F814W filter has an effective wavelength of 
A ~ 8090A and also probes the non-ionizing UV contin- 
uum at z ~ 3.09. The footprint of our NB3640 image 
contains 109 color-selected LBG candidates, of which 28 
(including 2 QSOs) have previous spectroscopically con- 
firmed redshifts z > 3.055, as well as 41 LBG candi- 
dates without previous spectroscopic confirmation. The 
footprint also conta ins 110 LAE candidates identified by 
INestor et all (|201lD . 

LAE candidates were selected using the B, V, and 
NB4980 data. The NB4980 filter covers redshifted HI 
Lya A1216A for galaxies with 3.05 < z < 3.12. We cre- 
ated a linear combination of the broad B and V images 
after scaling to a common photometric zero-point, such 
that our so-called "BV" image has an effective wave- 
length of A ~ 4980A. The BV - NB4980 color is there- 
fore a measure of Lya equivalent width for galaxies with 
3.05 < z < 312. T he LAEs in the sample described by 
INestor et all (j2011l) were required to have NB4980 < 26 
and BV — NB4980 > 0.7, which corresponds to a Lya 
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rest-frame equivalent width (REW) of > 20A. We also 
subtracted the BV image from the NB4980 image to cre- 
ate a so-called "LyA" image, which was used to define 
the centroids of Lya emission from the LAE candidates. 

2.2. Spectroscopic Observations 

We performed spectroscopy in the SSA22a field us- 
ing both shallow and deep observations. The purpose of 
the shallow observations was to obtain spectroscopic red- 
shifts for LBG and LAE photometric candidates, while 
the deep observations were intended to provide detailed 
information for objects with NB3640 detections. Multi- 
object spectroscopy was performed using the blue arm 
of th e LRIS dichroic spe ctrograph on Keck I (O ke et al.1 
119951: ISteidel et all 120041 ) . Nine shallow slit masks were 
observed over the course of five observing runs in 2009 
June, 2009 September, 2010 July, 2010 August, and 2011 
May. A single deep mask was observed in 2010 August. 
Typical exposure times for the shallow masks were 5400 
seconds (ranging from 3600 to 6270 seconds), while the 
deep mask was observed for 23700 seconds. For most 
runs, the conditions were photometric, with seeing rang- 
ing from 0'.'5 to 0'.'7. In 2011 May, when additional data 
were collected for one of the shallow masks, conditions 
featured variable cloudiness and seeing ranging from 0'.'8 
to l'.'O. 

The primary targets for shallow masks were LBG 
and LAE photometric candidates without previously- 
determined spectroscopic redshifts. LBGs with spectro- 
scopic redshifts were added as filler targets on the mask. 
A total of 50 LBGs and 114 LAEs were targeted on the 
shallow masks. Slits were centered on the detections in 
the NB4980 (i.e., Lya) for LAEs, and in the i?-band 
(i.e., rest- frame UV continuum) for LBGs. Seven out of 
nine shallow masks were observed using a 300 line mm -1 
grism blazed at 5000 A, along with the "d680" dichroic 
beam splitter, sending light with wavelengths bluer than 
~ 6800A to the blue arm of LRIS. The spectral resolu- 
tion for these masks was R = 530. The two remaining 
shallow masks were observed at higher spectral resolution 
using a 600 line mm" 1 grism blazed at 4000 A, one with 
the "d560" dichroic (splitting the incoming light beam at 
~ 5600 A), and the other with the "d680" dichroic. The 
spectral resolution for these two masks was R — 1200. 

The primary targets for the deep mask were LBGs and 
LAEs with NB3640 detections. A sample of 4 LBGs and 
13 LAEs with NB3640 detections were targeted on the 
deep mask (along with one LBG and 4 LAEs lacking 
NB3640 detections, but added as filler). Slits for objects 
with NB3640 detections were centered on the detections 
in the NB3640 image (i.e., LyC for objects at z > 3.055), 
while the NB4980 and _R-band detections were used, re- 
spectively, for centering the slits for the filler LAEs and 
LBG without NB3640 detections. The deep mask was 
observed using a 400 line mm -1 grism blazed at 3400 A, 
along with the "d680" dichroic beam splitter, sending 
light with wavelengths bluer than ~ 6800A to the blue 
arm of LRIS. The spectral resolution for this mask was 
R = 700. 

The data were primarily reduced using IRAF tasks, 
with scripts designed for cutting up the multi-object slit- 
mask images into individual slitlets, flat-fielding using 
spectra of the twilight sky, rejecting cosmic rays, sub- 



tracting the sky background, averaging individual expo- 
sures into final stacked two-dimensional spectra, extract- 
ing to one dimension, wavelength and flux calibrating, 
and shifting into the vacu um frame. T h ese pr ocedures 
are described in detail in [Stcidcl et al. (2003). There 
were a couple of notable differences in the data reduction 
procedures used for this sample, relative to the typical 
LBG reduction strategy. First, we used a custom IDL 
script (N. Reddy 2010, private communication) to rectify 
the curved slitlets before performing any of the standard 
IRAF reduction tasks. Since the majority of our targets 
are LAEs with negligible continuum and a single emis- 
sion line, we required slit rectification in order to extract 
a spectrum over a broad wavelength range at the loca- 
tion of the object indicated by the isolated bright Lya 
feature. Likewise, for deep-mask spectra, the faintness of 
the continuum level in the LyC region precluded a robust 
trace without rect ification. We also fol lowed the proce- 
dures outlined in iShaplev et al.l (|2006l ) for background 
subtraction of deep-mask spectra. Accordingly, to avoid 
potential over-subtraction of the background, the object 
continuum location was excluded from the estimate of 
the background fit at each dispersion point. We used 
the maximum possible number of pixels to fit the sky 
emission for each object. In practice, the widths of the 
sky regions on either side of the continuum location de- 
pended on the length of each slitlet and the position of 
the object along the slit. 

For LAEs, redshifts were calculated from the observed 
centroid of the Lya emission feature (rest- frame A Lya — 
1215.67 A). For LBGs, emission redshifts were estimated 
from the observed centroid of Lya, and absorption red- 
shifts from the centroids of interstellar metal absorption 
features when present. 

Finally, as described in Siana et al. (in preparation), 
near-IR spectra of several objects in the s ample were ob- 
taine d in 2011 August with NIRSPEC (jMcLean et al.l 
1998) on Keck II. These data were collected for a sep- 
arate project; we refer to these data here as, in three 
cases, the near-IR spectra reveal emission features rele- 
vant to the interpretation of the NB3640 detections (see 
Sections |4~1~21 and IIX3")) . 

3. THE UPDATED LBG AND LAE SAMPLES 

The INestor et all ()2011[ ) sample contained 26 LBGs 
with spectroscopic redshifts z > 3.055. Our new data 
include spectra of 11 LBGs and 1 QSO with previously 
determined redshifts, and 41 candidate LBG sources. For 
each of the 12 re-observed galaxies, the redshifts deter- 
mined from our new data agree with the previous deter- 
minations within Az = 0.012. The spectrum of one of 
the newly observed LBGs candidates reveals it to be a 
star. Of the 40 remaining galaxies, we were able to de- 
termine redshifts for 26 objects by identifying Lya emis- 
sion and/or interstellar metal absorption features. All 26 
galaxies have z > 2.45, 16 of which have z > 3.055. Thus, 
our updated sample of LBGs with redshifts z > 3.055 
contains 42 galaxies (26 previously known and 16 new). 
For this sample, the NB3640 filter probes the redshifted 
LyC spectral region. Coordinates, redshifts, and i?-band 
photometry are listed in Table [T] for the 42 z > 3.055 
sources, and in Table [2] for the 10 z < 3.055 newly con- 
firmed LBGs. 
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iNestor et all (|2011l ) identified 110 LAE candidates 
having BV - NB4980 > 0.7 and NB4980 < 26. Our new 
data include spectra of 96 of the 110 LAE candidates. 
Four of the LAEs for which we did not acquire new spec- 
tra are also LBGs with previously determined spectro- 
scopic redshifts: C4, C28, M28 and MD23. The other 10 
LAE candidates for which we did not acquire a spectrum 
are relatively faint in NB4980 (25.25 < NB4980 < 26) 
but are randomly distributed in BV — NB4980 color. 
Thus, this incompleteness should not bias our results. 
We detect an emission line in the expected spectral re- 
gion, « 4935A — 5015A, in the spectra of 88 of the can- 
didates. One of the 88 is D3, which is also in our LBG 
sample with a previous spectroscopic redshift, and an- 
other two, C9 and M13, were LBG candidates that are 
now spectroscopically confirmed members of our LBG 
sample. Thus, there are seven objects that appear in 
both our LBG and LAE samples. 

In principle, some of the emission lines detected in 
our spectra could be [OII]A3727, or in some cases H/3 
or [OIII]A5007, from very faint lower- redshift systems. 
However, 76 of the 88 (86%) of the sources with a line 
detection have spectra with sufficient spectral coverage 
that, if the detected line was [Oil], H/3, or [OIII] at lower 
redshift, at least one of the other rest-frame optical fea- 
tures should have been detected as well. In no cases 
do we detect such corresponding lines. Furthermore, in 
a similar LAE survey at z = 3.1 in th e Exte nded Chan- 
dra Deep Field-South. iGronwall et all (|2007l ) argue that, 
given the relatively small volume covered by the narrow- 
band filter at z — 0.34 relative to z = 3.1 and the rarity of 
[Oil] emitting galaxies with REW above their threshold 
(REW> 60A at z = 0.34), their level of contamination by 
low- 2; [Oil] emitters is ne gligible. As our filter i s ~ 60% 
broader than that used bv lGronwall et all (|2007l ). we rec- 
ognize a slight possibility of misidentification of a small 
number of emission lines. However, we continue with the 
assumption that all 88 of the lines that we determined 
to be HI Lya are correctly identified as such. One of 
these spectroscopically confirmed objects, LAE034, has 
z = 3.044. At this redshift some non-ionizing UV flux 
will contribute to the NB3640 detection, and therefore 
we excluded it from our sample. 

Of the eight candidates for which we did not detect any 
emission line in the « 4975 ± 40A spectral region, seven 
are relatively faint in NB4980, and the other is very dif- 
fuse in the LyA image with an extent significantly larger 
than the slit width. Thus, these objects may also be 
Lya emitting galaxies at z > 3.055 that have insufficient 
line fluxes and/or sufficient slit losses such that the Lya 
emission line is not detectable in our data. In order to 
quantify the expected detection significance levels of Lya 
emission lines in our spectra, we predicted the Lya line 
fluxes for each LAE candidate using the measured BV 
and NB4980 magnitudes. We then measured the noise 
properties of the calibrated spectra in the wavelength in- 
terval corresponding to the width of the NB4980 filter 
in order to assess the minimum detectable line fluxes as- 
suming unresolved lines. In practice, this method will 
underestimate the detection limit in some of our spectra, 
as we do not quantitatively account for slit losses or the 
possibility of resolved line profiles. 

Figure Q] shows the predicted detection significance 
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Fig. 1. — The predicted detection significance level (SL) for 
emission lines in the spectra of our LAE candidates. The SL 
values are calculated using the photometrically-estimated line 
fluxes and the noise properties of the spectra. They do not 
account for slit losses or the possibility of resolved lines, and 
are thus upper limits. The predicted SL values for spectra 
with and without a detected emission line are represented 
by the black dots and red stars, respectively. Only one of 
our non-detections, LAE020, has a relatively high predicted 
SL value. LAE020 appears very diffuse in the LyA image, 
however, and thus its spectrum is expected to suffer from 
large slit losses. 

level (SL), i.e., the ratio of the photometrically-predicted 
line flux to the approximated line flux uncertainty, for 
each LAE candidate observed. For objects with mul- 
tiple spectra having different resolutions, we show the 
smallest SL value if Lya is detected, and the largest SL 
value if it is not. We do not consider spectra in which 
the slit did not cover the position of the LyA flux, as 
is the case for some of the deep-mask spectra for which 
the slit was positioned over the offset NB3640 flux. All 
of the detections are consistent with having predicted 
SL > 3, with a median value of SL = 13. Seven of the 
eight objects without a detected emission line have, rela- 
tive to those with detections, small upper-limits to their 
predicted detection significance, SL < 6, with a median 
value of SL = 4; the other object is the diffuse system 
(LAE020) mentioned above. Thus, we can not rule the 
non-detections out as being z ~ 3.09 LAEs. Addition- 
ally, considering the high success rate (> 91%) of our 
photometric selection at identifying z > 3.055 LAEs, it 
is likely that most of the 10 candidates for which we have 
no spectroscopic data are also at z > 3.055. Nonetheless, 
we conservatively removed the 18 candidate LAEs with- 
out spectroscopic confirmation from the statistical sam- 
ple discussed in this work, leaving a spectroscopic LAE 
sample of 91 galaxies spanning 3.057 < z < 3.108. We 
summarize our photometric and spectroscopic results for 
t he LAE samp le in Ta ble |31 

INestor et all (|2011D also identified an additional 20 
"faint sample" LAE candidates with 26 < NB4980 < 
26.5 and BV - NB4980 > 1.2, corresponding to LAE 
IDs 111 — 130. We obtained spectra for 18 of these faint 
LAE candidates. We were unable to determine redshifts 
for two such systems (again, in data with low SL val- 
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Fig. 2. — LBGs with z > 3.055 detected in NB3640. Images are centered on the 7?-band centroid and span 7"x7". The 
orientation is such that north is up and east is to the left. At these redshifts, the NB3640 filter (image shown here after 
smoothing by a Gaussian kernel with FWHM=0'.'5) samples the rest-frame LyC, while the R and F814W filter s sample the 
rest-fr ame non-ionizing UV continuum. The green contours indicate NB3640 flux levels. Left: LBGs discussed in Nestor et al.l 
(|20Tlh . Right: Newly confirmed z > 3.055 LBGs. 



ues). All of the other 16 systems were determined to 
have 3.069 < z < 3.111. The faint sample LAEs have 
larger photometric uncertainties than the main sample 
LAEs and w ere selected with sli ghtly different photomet- 
ric criterea (jNestor et all 1201 If ). Thus, they are not in- 
cluded in the statistical analyses presented in this work. 
We include them in Table [3] for completeness. 

4. SYSTEMS WITH NB3640 DETECTIONS 

With our updated spectroscopic sample of LBGs and 
LAEs in place, we now consider the set of objects with 
detections in our NB3640 imag e . Of the 26 z > 3.055 
LBGs discussed in Nes tor et all (|2011h . six are detected 
in NB3640. Four of our 16 newly confirmed high-redshift 
LBGs are also detected in NB3640: C14, M2, M5, and 
M29. Their NB3640-, R-, and HST F814W images are 
shown in Figure [21 with the contours corresponding to 
28.81, 28.06 and 27.62 mag arcsec" 2 in the NB3640 im- 
age. All four of the newly confirmed LBGs with NB3640 
detections appear clumpy in the H ST / ACS images, with 
multiple discrete regions of non-ionizing continuum. The 
NB3640 flux appears to cover all of the clumps in C14, 
M2 and M5. In M29, the NB3640 flux appears to be 



associated with only the more compact clump, which is 
~ 0'.'75 to the northeast of a more extended region of 
non-ionizing UV flux. With the addition of these four 
new NB3640 detections, we now have 10 NB3640 detec- 
tions associated with 42 z > 3.055 LBGs. Their NB3640 
magnitudes, the spatial offset between the i?-band and 
NB3640 flux centroids, and the inferred non-ionizing to 
ionizing UV flux-density ratios are presented in Table [TJ 
In our sample of 91 spectroscopically confirmed LAEs, 
20 have N B3640 detections. A n additional 6 LAE candi- 
dates from lNestor et all (|2011D have NB3640 detections: 
one was not targeted in our spectroscopy, and we were 
unable to confirm the redshifts of five other candidates. 
Thus, these six objects are not discussed here. For the 
20 LAEs with associated NB3640 detections, we report 
in Table [3] the NB3640 magnitudes, spatial offsets be- 
tween both i?-band and LyA and NB3640 flux centroids, 
and the inferred non-ionizing to ionizing UV flux-density 
ratios. 

Below, we discuss in detail the combined spectroscopic 
and imaging data sets for the individual objects in our 
updated samples that were considered LyC detections 
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Fig. 3. — The emission structure of MD46. Left: The spa- 
tial extent of Lya emission compared to the rest-frame non- 
ionizing UV continuum emission and putative LyC radiation 
(green contours). Panel (a) shows the -ffS"T/ACS-F814W im- 
age (spatial resolution ~ O'.'l) while panel (b) is our "LyA" 
(seeing ~ 0'.'8) image. The slit position for our deep-mask 
spectrum is indicated in both images. Panel (c) displays the 
spatial (i.e., along the slit) extent of the Lya line, averaged 
over (rest-frame) ~ 4A of our deep-mask spectrum and reg- 
istered to the images. The dotted curve is a Gaussian fit to 
the profile, excluding the excess at positive offset. Right: The 
relative spatial extents of LyC and Lya emission. Panel (d) 
shows the two-dimensional spectrum in the redshifted Lya 
region. The two extraction apertures applied in Figure [4] are 
indicated with cyan boxes. In panel (e), we show the profile 
for the region of the two-dimensional spectrum averaged in 
the spectral direction from rest-frame A ~ 875A to A ~ 9I0A 
(green), as well as that of the Lya region, re-scaled for ease 
of comparison (black). 

in INestor et al.l (| 2 1 If) , as well as the newly determined 
z > 3.055 LBGs. The data are shown in Figures [3TfT5l 
ifST/ACS-F814W images are shown when available, 
BV images are shown otherwise, and all of the images 
span 8" per side. Slit positions are indicated by blue 
(deep mask) and/or black (shallow mask) boxes. Green 
contours represent NB3640 flux levels (see above). The 
LAE images also contain red contours, representing LyA 
(see Section 12. ip flux levels^ Two-dimensional spec- 
tra, when shown, are registered in the spatial direction 
to match the accompanying imaging and span ~ 125A 
in the spectral direction. Unless otherwise noted, one- 

3 Due to the large range of LyA fluxes in our LAE sample, 
the flux levels represented by the red contours vary from image to 
image. 
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Fig. 4. — The extracted deep-mask spectra of MD46. We 
identify the emission feature at A = 4973A as Lya at red- 
shift z = 3.091 (vertical dashed line). The corresponding Ly- 
man break is marked with a vertical dotted line. The black 
spectrum, which has been boxcar smoothed (5 pixels), is the 
extraction containing only the northwestern clump (see Fig- 
ure^, which is spatially consistent with the detected NB3640 
flux. The red spectrum represents the extraction that ex- 
cludes the northwestern clump. The red spectrum has been 
scaled such that the average flux levels redward of Lya are 
approximately equal. Also shown in green, below the Lyman 
limit, is the NB3640 filter transmission curve, 
dimensional spectra are shown without smoothing or bin- 
ning. When multiple one-dimensional spectra are shown, 
the additional spectra are offset for clarity and the non- 
offset spectrum corresponds to the horizontal slit in the 
imaging (by default, the deep-mask slit when available). 

4.1. LBGs 

In this section, we discuss the individual LBGs hav- 
ing NB3640 detections for which we have new data. We 
begin with the four systems which we retain as possible 
LyC-leaking galaxies. We next discuss the three LBGs 
for which we find evidence for the presence of a fore- 
ground interloper in our new data. We then discuss the 
special case of aug96M16, and conclude the section with 
a summary of our LBG NB3640 detection sample. 

4.1.1. LBG Lyman- continuum Candidates 
MD46 

The LBG MD46 (z = 3.091) has a complex multi- 
component morphology, which can be clearly seen in Fig- 
ure 02 Panel (a) shows the FST/ACS-F814W image, in 
which the flux appears to originate from at least three 
distinct clumps. Panel (b) shows the ground-based LyA 
image. The clump to the east, which partly enters the 
slit, is spatially coincident with negative flux in the LyA 
image (corresponding to Lya absorption at z ~ 3.09). 
The bulk of the Lya emission appears to emanate from 
in and around the central clump, while the bulk of the 
NB3640 flux avoids this central clump and is instead spa- 
tially coincident with the clump to the northwest, 
which also falls on the slit. Thus, it is important to de- 
termine if the northwestern clump is also at z ~ 3.09. 
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Fig. 5.— HST/ ACS F814W images (8" per side) and spectra of LBGs with NB3640 detections retained as possible LyC 
emitters. The positions and sizes of the slits are indicated with boxes, and green contours indicate NB3640 flux levels. When 
present, vertical dashed and dotted lines in the spectra indicate the position of the redshifted Lya line and Lyman break, 
respectively. Top: Image and shallow-mask spectrum C14. The other shallow-mask spectrum of C14 (corresponding to the 
diagonal slit) has a very low signal to noise ratio and is not shown. Middle: The F814W image and combined shallow-mask 
spectrum of M2. Bottom: Image and shallow-mask spectrum of M29. The locations of several common interstellar absorption 
features are marked with vertical dashed (detected) and dash-dotted (non-detected) lines. 



The region of the two-dimensional spectrum corre- 
sponding to redshifted Lya emission is shown in panel 
(d). The Lya emission line appears asymmetric, which 
can be seen more clearly in panel (c), where we have av- 
eraged over r; 4A (rest-frame) in the spectral direction, 
centered on the peak of the Lya emission line, and fitted 
a Gaussian profile (dotted curve) constrained only by the 
data south of the spatial zero-point indicated by the fig- 
ure axis. This fit highlights the asymmetric northwestern 
extension, which is spatially coincident with the north- 
western clump seen in the HST imaging and the NB3640 
flux in our Keck imaging. The spatial coincidence of the 
Lya and NB3640 emitting region implies that the north- 
western clump is also very likely to be emitting Lya at 



z ~ 3.091, which in turn indicates that the NB3640 flux 
is indeed escaping LyC emission. 

Panel (e) shows the spatial profile for the region of 
the two-dimensional spectrum just below the redshifted 
Lyman-limit (green histogram), averaged in the spectral 
direction over A ~ 880A — 910A in the rest-frame, as well 
as that of the Lya region (black histogram, re-scaled by 
factor of 26 for ease of comparison). Again, the NB3640 
flux is spatially aligned with the extended Lya, suggest- 
ing that it is actually LyC flux escaping from a galaxy at 
z = 3.09. We extracted one-dimensional spectra in two 
locations in the two-dimensional spectrum. The extrac- 
tion apertures are shown in panel (d). Figure [4] displays 
the resulting one-dimensional spectra. The black spec- 
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Fig. 6. — As in Figure [5] but for systems with evidence of contamination (top and middle) or with unconfirmed redshift 
(bottom). Top: Image and deep-mask spectrum of MD32. Middle: Image and shallow- (offset, upper) and deep-mask (lower) 
spectra of M5. Bottom: Image of aug96M16 and the deep-mask spectrum. We were unable to determine the redshift of the 
region most closely associated with the NB3640 emission. 



trum, which has been smoothed with a 5-pixel boxcar, 
represents the northwestern extraction corresponding to 
the extended Lya emission and the detected NB3640 
flux. The red spectrum, which has been scaled down 
by a factor of 5.6 such that the average flux level red- 
ward of Lya emission is equal to that of the black spec- 
trum, represents the southern extraction, corresponding 
to the central clump, which is undetected in our NB3640 
imaging. The spectrum of the central clump appears 
to exhibit a strong break at the Lyman limit, the posi- 
tion of which is indicated by the vertical dashed line. In 
contrast, the northwestern clump appears to lack a sharp 
break at the Lyman limit. We also note that several Lya 
forest absorption lines appear to be present in both spec- 
tra (e.g., at A ~ 4720A, 4860A, 4925A, etc.) suggestive 
of the continua originating at similar emission redshifts. 
However, the poor signal to noise ratios prevent quanti- 



tative analysis of any perceived correlation. 

We conclude that the z = 3.091 LBG MD46 contains 
three distinct (rest-frame) UV-bright regions. The south- 
eastern region is a Lya absorber as seen from our viewing 
perspective, while the central and northwestern regions 
exhibit Lya in emission. Of these latter regions, the 
fainter region (in both continuum and line emission) is 
more compact and can be seen in ionizing continuum 
along our sightline. Thus we retain MD46 in our LyC- 
emitting LBG sample. 

C14 

We obtained spectra of the LBG candidate C14 on two 
of the shallow masks. Each spectrum exhibits an emis- 
sion line at A ~ 5130A, which we attribute to Lya at 
z = 3.220. C14 is detected in our NB3640 image (Fig- 
ure [2]). In the H ST /ACS image, shown in the top panel 
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Fig. 7.— The emission structures of (from top down) LAE010, LAE016, and LAE018. The left panels show the HST/ACS- 
F814W images. The boxes indicate the positions of all shallow- (blue) and deep- (black) mask slits. Green and red contours 
represent flux levels in the NB3640 and LyA images, respectively. The right panels display the one-dimensional extracted 
spectra. The vertical dashed/dotted lines indicate the position of the redshifted Lya/Lyman limits, respectively. For LAE018, 
we show both the shallow- (offset, upper) and deep-mask (lower) spectra. 



of Figure O C14 breaks into a brighter and a fainter 
clump, with the NB3640 emission appearing to span the 
two clumps. As no other emission features were robustly 
detected in either spectrum, we retain C14 as a possi- 
ble LyC-emitting galaxy. We note, however, that the 
spectrum corresponding to the slit position that is bet- 
ter aligned with both clumps (not shown) is of relatively 
poor quality, and therefore our limits on the absence of 
other emission features from the fainter clump are rela- 
tively weak. 

M2 

We obtained spectra of the LBG candidate M2 on two 
of the shallow masks. Both spectra exhibit an emis- 
sion line at A ~ 5334A, which we attribute to Lya at 
z = 3.388. M2 is detected in our NB3640 image (Fig- 



ure H]). In the HST/ ACS image, shown in the middle 
panel of Figure M2 breaks into a brighter and a fainter 
clump, with the NB3640 emission appearing to span the 
two clumps. No other emission features were robustly 
detected in either spectrum. Thus, we retain M2 as a 
possible LyC-emitting galaxy. 

M29 

We obtained a spectrum of the LBG candidate M29 
on one of the shallow masks. The spectrum does not 
exhibit emission lines. However, we identify interstel- 
lar absorption lines from OI, SilV, Sill, and CIV corre- 
sponding to a redshift of z = 3.228. M29 is detected in 
our NB3640 image (Figure [5]) . Two clumps of emission 
are seen in the H ST/ ACS image, shown in the bottom 
panel of Figure [5l The NB3640 emission only covers the 
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Fig. 8.— Same as Figure but for LAE021, LAE028, and LAE038. Our BV image is shown for LAE021 as it does not 
have H ST/ ACS-F814W imaging. For LAE028 and LAE038, we show both the shallow- (offset, upper) and deep-mask (lower) 
spectra. 



more compact clump. It is possible that both clumps 
are at z = 3.228, with the NB3640 emission emanating 
from a compact star-forming clump. Alternatively, the 
more diffuse source may be at high redshift while the 
compact source associated with the NB3640 emission is 
in the foreground. Lacking any direct evidence to favor 
either scenario, we retain M29 as a possible LyC-emitting 
galaxy and account for the possibility of contamination 
with our Monte Carlo simulation (Section [5]). 

4.L2. LBGs With Evidence For Contamination 

MD32 

The HST / ACS image and our deep-mask spectrum of 
MD32 (z — 3.102) are shown in the top panel of Fig- 
ure [SJ Our spectrum exhibits flux below the redshifted 
Lyman limit, and contains no evidence for foreground 
contaminants. However, in the NIRSPEC observations 



the slit was aligned along the direction of elongation of 
MD32. The resulting two-dimensional spectrum exhibits 
an additional line within ~0'/7 of MD32 which, if [OIII] 
A5007, would correspond to a redshift of z = 2.88. Due 
to this evidence of contamination, we remove MD32 from 
our list of possible LyC-emitting galaxies. 

M5 

We obtained spectra of the LBG candidate M5 with 
both the deep mask and one of the shallow masks. Each 
spectrum exhibits an emission line at A ~ 5260A, which 
we attribute to Lya at z = 3.327. M5 is detected in our 
NB3640 image (Figure [2]). Two clumps of emission are 
observed in the HST / ACS image, shown in the middle 
panel of Figure [B] The NB3640 emission appears to span 
the two clumps. However, our deep-mask spectrum of 
M5, also shown in Figure [6] along with the shallow-mask 
spectrum, reveals the presence of an additional, bluer, 
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Fig. 9.— Same as Figure but for LAE039, LAE041, and LAE046. BV images are shown as LAE039, LAE041 and LAE046 
do not have ST/ACS-F814W imaging. For LAE041, we show both the shallow- (offset, upper) and deep-mask (lower) spectra. 



emission line at A ~ 4040A, indicating the presence of 
a lower-redshift interloper. It is noteworthy, however, 
that the detected NB3640 flux covers both clumps of 
continuum emission. Thus, if the interloper is associated 
with only one of the clumps, it may be that M5 is in- 
deed being detected in LyC emission. Nonetheless, we 
conservatively remove M5 from our list of possible LyC 
detections. 

C49 

Although we did not obtain a new LRIS spectrum of 
C49 (z = 3.163), it was observed with NIRSPEC. As 
shown in Figure [2j the HST / ACS image reveals that 
C49 comprises two distinct clumps. The NIRSPEC slit 
was positioned to capture the spectra of both of these 
clumps. While the spectrum of the southern clump con- 
firms the redshift z = 3.16, that of the northern clump, 
which is spatially closer to the NB3640 detection, has an 



additional emission line which, if [OIII] A5007, would cor- 
respond to a redshift of z = 2.97. We therefore remove 
C49 from our list of possible LyC-emitting galaxies. 

4.1.3. Unconfirmed LBG redshift 

aug96M16 

The HST /ACS image of aug96M16 is shown in the 
bottom panel of Figure M As with MD46, the HST 
imaging reveals the emission from aug96M16 to comprise 
several distinct clumps. The detected NB3640 flux is co- 
incident with the central region of emission that spans 
the l'/2 width of our LRIS slit. Only the western-most 
clump (w 1'.'3 west of the slit edge) was resolved as a dis- 
tinct clump, separate from LBG candidate aug96M16, in 
the ground-based data originally used to select LBG can- 
didat es based on U — G and G — 1Z colors in the SSA22a 
field ( Steidel et al. 2003). This western-most clump has 
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Fig. 10.— Same as Figure [7] but for LAE048, LAE051, and LAE053. For LAE053, we show both the shallow- (offset, upper) 
and deep-mask (lower) spectra. Lya appears in absorption in the deep-mask spectrum of LAE053, as the slit was centered on 
the NB3640 detection which is offset by ~ 0'.'9 from the LyA emission. 



a previously determined spectroscopic redshift z = 3.285, 
which was al s o (err oneously) attributed to aug96M16 in 
I Nest or et al.l ()2011l ). The slit position in our LRIS deep 
mask covered neither the western-most nor eastern-most 
clump (« 0'.'75 east of the slit edge). The slit position of 
our NI RSPE C spectrum (Siana et al., in preparation; see 
Section l2~2j) of aug96M16 was aligned to cover all of the 
emission clumps, however. Using the NIRSPEC data, we 
determine a redshift z — 3.09 for the eastern-most clump 
and confirm the redshift of the western-most clump as 
z = 3.29. The large difference in redshift between the the 
eastern-most and western-most clumps indicates that, al- 
though both are at z > 3.055, they are physically unre- 
lated. Our deep-mask spectrum of the central clumpy 
region, which is coincident with the NB3640 detection, 
is also shown in Figure |5] We are unable to determine 
a redshift for this region from either our deep mask or 
NIRSPEC spectra. Therefore, although it was photo- 



metrically selected as an LBG candidate and is within 
~ 1" on the sky to a spectroscopically-confirmed galaxy 
at z = 3.09, we conservatively remove aug96M16 from 
our list of possible LyC-emitting galaxies as well as from 
the parent z > 3.055 LBG sample. 

4.1.4. Summary of LBGs 

As discussed above, our sample began with 10 LBGs 
with possible LyC detections. Our new data includes 
spectra of seven of these 10 objects. In these spectro- 
scopic data we find evidence for the presence of a fore- 
ground object in close proximity on the sky to three of 
these LBGs - MD32, C49, and M5 - and are unable 
to confirm the redshift of the source associated with the 
NB3640 flux in a fourth - aug96M16. Of the seven LBGs 
having NB3640 detections with new spectroscopic data, 
we retain as possible LyC emission NB3640 detections in 
four: MD46, C14, M2, and M29, as well as two detections 
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Fig. 11— Same as Figure [7] but for LAE064, LAE069, and LAE074. Our BV image is shown for LAE069 as it does not 
have H 5T/ACS-F814W imaging. For LAE074, we show both the shallow- (offset, upper) and deep- mask (lower) spectra. Lya 
emission is not detected in the deep-mask spectrum. As the slit was centered on the NB3640 detection the bulk of the LyA 
emission fell outside of the slit. 



for which we do not have new data: D17 and C16. Thus, 
our LBG sample now contains of a total of six putative 
LyC-leaking galaxies from a parent sample of 41 LBGs 
(i.e., excluding aug96M16; see Section [4. 1 .3[) . Notably, 
we are able to study the spatial distribution of the Lya 
and NB3640 emission in detail for MD46, and find com- 
pelling evidence that the NB3640 emission is escaping 
LyC flux. 

It is possible that some of the 32 LBGs with no NB3640 
detections have levels of escaping LyC flux that are below 
our detection limit of tonb3640 ~ 27.3. To investigate 
this possibility, we stacked cut outs of the NB3640 image 
centered on the locations of these 32 LBGs. We detect 
no flux in the stacked image down to a 3cr limit on the 
average magnitude of mNB3640 = 27.96. 

4.2. LAEs 



As above, in this section we discuss the individual 
LAEs with NB3640 detections. We begin with an 
overview of the 17 systems which we retain as possible 
LyC-leaking galaxies. We next discuss the three LAEs 
for which we find evidence for the presence of a fore- 
ground interloper. We then discuss the six LAEs with 
NB3640 detections for which we were unable to confirm 
redshifts and are thus not included in our sample, and 
conclude the section with a summary of our LAE NB3640 
detections. 

4.2.1. LAE Lyman- continuum Candidates 

As indicated in Table EJ we spectroscopically confirm 
redshifts of 3.070 < z < 3.108 with no evidence for con- 
tamination for 17 LAEs with NB3640 detections. The 
images and spectra for these LAEs are shown in Fig- 
ures [7] — [I2j We briefly discuss each individual field in 
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Fig. 12.— Same as Figure but for LAE081 and LAE083. For LAE081, we show both the shallow- (offset, upper) and 
deep-mask (lower) spectra and BV image as it does not have H ST / ACS-F814W imaging. 



the appendix. Of particular note is LAE053, for which 
the LyA and NB3640 emission are offset by ~ 0'/9. The 
F814W image reveals that each of the LyA and NB3640 
detections is coincident with one of two distinct emission 
clumps. Our shallow- and deep-mask slits cover only 
the clump associated with the LyA and NB3640 fluxes, 
respectively. In the shallow-mask spectrum, we detect 
Lya in emission at z = 3.090. In the deep-mask spec- 
trum, however, we detect an absorption feature at the 
same wavelength as the emission line in the shallow-mask 
spectrum. We interpret this feature as Lya absorption at 
z = 3.090, which is compelling evidence that the clump 
associated with the NB3640 emission is also at z = 3.090 
and, in turn, that the NB3640 flux is indeed LyC emis- 
sion. In contrast, however, we caution that for four of the 
candidates shown in Figures I7HT21 (LAE021, LAE046, 
LAE048 and LAE069), the position of the NB3640 flux 
was not covered by any of the slit spectra. Thus we are 
unable to search directly for evidence of foreground con- 
tamination being responsible for the detected NB3640 
flux. Additionally, the F814W images of LAE046 (Fig- 
ure© and LAE048 (Figure [TO]) reveal that the NB3640 
and LyA fluxes, respectively, are coincident with sources 
of non-ionizing continuum that are significantly (> 1'.'5) 
spatially offset from each other, suggesting the sources 
are unrelated. In the absence of direct spectroscopic evi- 
dence confirming contamination for these individual sys- 
tems, we ran a Monte Carlo simulation (described in Sec- 
tion [SJ to statistically address the possibility of contami- 
nation in the entire sample of 17 detections. In this sim- 
ulation, NB3640 detections with large offsets are likely to 
be rejected as foreground interlopers. Thus, while indi- 
vidually the conclusions that LAEs such as LAE046 and 



LAE048 are LyC-leaking galaxies should be taken with 
caution, our statistical results are robust to the possibil- 
ity of foreground contamination. 

4.2.2. LAEs With Evidence For Contamination 

We find evidence for the presence of a foreground 
galaxy in the spectra of three of our LAEs, which we 
prese nt below. Additiona lly, one of the LAE candidates 
from iNestor et all (1201 ID . LAE034, is at a low enough 
rcdshift that it experiences contamination of its NB3640 
flux from non-ionizing UV radiation. We have thus re- 
moved LAE034 from our list of possible LyC-emitting 
galaxies as well as from the parent z > 3.055 LAE sam- 
ple. 

LAE003 

LAE003 has the brightest NB3640 detection 
(ffiNB364n = 24.74 ) of any of our LBGs or LAEs in 
INestor et all (|2011ft . The top-left panel of Figure Q3] 
shows the F814W image indicating the relative positions 
of the non-ionizing UV flux, the deep and shallow-mask 
slits, the Lya flux and the NB3640 flux. The middle 
panel shows the region of the two-dimensional spectrum 
corresponding to the expected location of the redshifted 
Lya line. We identify the line detected at A = 4980 as 
Lya, indicating that LAE003 is at z = 3.097. 

However, as can be seen in the lower panel, which 
shows a bluer region of the two-dimensional spectrum, 
there is another emission line spatially consistent with 
the LAE. If this line is also Lya emission, it corresponds 
to a redshift of z = 1.763. The bluer line is very slightly 
offset to the south, consistent with the tail-like structure 
seen in the H ST / ACS image. It is therefore possible that 
the z = 3.097 LAE003 is opaque below the Lyman limit 
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Fig. 13.— The emission features of LAE003. Left: The top panel shows the ffi?T/ACS-F814W image. The middle panel shows 
the two-dimensional spectrum in the redshifted Lya region, while the bottom panel shows a redder region of the spectrum, 
centered at the wavelength (3359A) of a bluer emission line. Right: The extracted one-dimensional spectra of LAE003. Shown 
are both the shallow- (offset, upper) and deep-mask (lower) spectra. We identify the emission feature at A = 4984A as Lya at 
redshift z = 3.097. The corresponding Lyman break is marked with a vertical dotted line. The bluer emission line is clearly 
detected at A = 3359A in the deep-mask spectrum. The green curve below the Lyman limit indicates the shape of the NB3640 
filter. 
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Fig. 14. — The emission features of LAE019. Left: The top panel shows our BV image. The middle and bottom panels show 
the two-dimensional deep-mask spectrum in the regions of the redshifted Lya, and the bluer emission line, respectively. Right: 
The extracted one-dimensional spectra of LAE019. Shown are both the shallow- (offset, upper) and deep-mask (lower) spectra. 
We identify the emission feature at A = 4985A as Lya at redshift z = 3.101. The corresponding Lyman break is marked with a 
vertical dotted line. The bluer emission line can be seen at A = 3490A in the deep-mask spectrum. 
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LAE025 W 




Fig. 15. — The emission features of LAE025. The top panels 
show the FST/ACS-F814W image. The boxes indicate the 
positions of the slits in the shallow- (left) and deep- (right) 
masks. The inset in the top-left panel shows the central region 
without contours after slight (2.35-pixel FWHM) smoothing. 
The middle panels show the two-dimensional spectra in the 
redshifted Lya region for the corresponding slit spectra, with 
an emission line at ~ 4974A in the shallow-mask spectrum 
and, possibly, at ~ 4942A in the deep-mask spectrum. The 
bottom panels show the two-dimensional spectra in the vicin- 
ity of the emission line at ~ 4035A. 

and the NB3640 flux is non-ionizing continuum emerg- 
ing from this lower-redshift galaxy. It is noteworthy that 
the NB3640 flux is not centered on the portion of the 
F814W flux that we associate with the lower-redshift in- 
terloper, and extends to the north-west. However, due 
to the clear presence of an interloper we take the con- 
servative approach and remove LAE003 from our list of 
possible LyC-leaking galaxies. The one-dimensional ex- 
tracted spectra are shown in the right panel of Figu rc ll3l 
LAE 003 had previously been studied by llnoue et all 
(|2011l ) who determined a spectroscopic redshift of z = 
3.100. Their spectrum did not extend blueward enough 
to detect the emission line at A = 3359 A, however. As 
the NB3640 (NB359 in llnoue et aTll2011h flux is spatially 
coincident with the i?-band and LyA flux for LAE003, 
these authors argue that the probability of foreground 
contamination is small and therefore the NB359 flux is 
indeed LyC. They interpret the implied very high ratio 
of ionizing to non-ionizing UV flux density of four ob- 
jects, including LAE003, in terms of very young stellar 
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Fig. 16. — The extracted shallow- (offset, upper) and deep- 
mask (lower) spectra of LAE025. We identify the emission 
feature at A = 4974A in the shallow-mask spectrum as Lya 
at redshift z — 3.091. The corresponding Lyman break is 
marked with a vertical dotted line. The bluer emission line 
at A = 4035A can be seen in both spectra. The deep-mask 
spectrum appears to exhibit Lya emission at A = 4942 A, 
corresponding to z = 3.065. Also shown in green, below the 
Lyman limit, is the NB3640 filter transmission curve. 

populations with top-heavy initial mass functions. This 
interpretation is invalidated for LAE003 by the discovery 
of the low-redshift interloper in our deep- mask spectrum. 
The case of LAE003 highlights the importance of obtain- 
ing spectra extending as far to the blue as possible for 
properly interpreting possible LyC-leaking galaxies. 

LAE019 

We do not have HST imaging of LAE019. The top 
panel of FigurefMlshows our BV image with the deep and 
shallow-mask slits indicated. The Lya emission centroid 
is slightly offset to the north in both the imaging and two- 
dimensional deep-mask spectrum (middle panel), while 
the NB3640 flux centroid is slightly offset to the south 
in our imaging. The deep-mask spectrum of LAE019, 
shown in Figure [T4l also exhibits a bluer emission feature 
at A ~ 3490A. If this feature is Lya emission, it would 
indicate the presence of an interloper at z = 1.872. The 
bottom panel shows the two-dimensional spectrum in the 
A ~ 3490A region. This lower-redshift emission line is 
spatially consistent with the NB3640 flux, indicating that 
the NB3640 flux is likely due to the interloper. Therefore, 
we remove LAE019 from our list of possible LyC-leaking 
galaxies. 

LAE025 

We obtained spectra of LAE025 in both shallow and 
deep masks. The top panels of Figure [TS] show the 
HST / ACS image, with the corresponding slit positions 
(shallow left; deep right). The inset in the top- left panel 
shows the central region without contours, after smooth- 
ing by a Gaussian kernel with FWHM= 2.35 pixels (O'.'l). 
The bulk of the NB3640 flux is associated with the cen- 
tral clump. We detect flux in neither the F814W nor BV 
images at the location of the LyA flux centroid. The 
middle-left panel of Figure [15] shows the region of the 
two-dimensional shallow-mask spectrum corresponding 
to the expected location of the redshifted Lya line. We 
detect a strong emission line spatially coincident with the 
LyA flux, which we identify as Lya at z = 3.091. How- 
ever, the non- ionizing UV flux detected in the H ST / ACS 
image that is coincident with the NB3640 flux is just 
off of the slit. The middle- right panel shows the same 
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spectral region of the deep-mask spectrum, for which the 
slit was centered on the NB3640 detection. No obvious 
emission at A ~ 4974A is seen, as might be expected if 
the region coincident with the NB3640 flux was also at 
z = 3.091. Based on possible detections of Lya emis- 
sion at A ~ 4942A and perhaps corresponding Sill in- 
terstellar absorption, we tentatively assign a redshift of 
z = 3.065. This is a high enough redshift such that the 
NB3640 filter is still opaque to radiation longward of the 
Lyman limit. It should be noted, however, that the LAE 
at z — 3.091 for which we searched for corresponding 
NB3640 flux would be unrelated to the z = 3.065 object 
(the relative velocity between the two Lya lines being 
~ 1800 kms -1 ) associated with the NB3640 flux; their 
proximity on the sky would be coincidental. 

Further complicating this system is the presence of an 
additional, much bluer line at A ~ 4035A, which is de- 
tected in both spectra (Figure [T6f . It is spatially consis- 
tent with the faint, low surface-brightness flux detected 
in the F814W image to the north- northwest of the cen- 
tral clump, which falls in both slits. If identified as Lya, 
it implies a redshift of z — 2.319. 

With the data at hand it is difficult to either confirm or 
refute the claim that the detected NB3640 flux is ioniz- 
ing continuum. Nonetheless, considering (a) there is ev- 
idence for a nearby foreground object, and (b) our iden- 
tification of the redshift of the blob associated with the 
NB3640 emission is tentative, we remove LAE025 from 
our list of possible LyC detections. 

4.2.3. Unconfirmed LAE Candidates with NB3640 
Detections 

We were unable to confirm the rcdshifts of six of 
our LAE candidates with NB3640 detections. Of these 
six, we did not observe one, LAE084. The other five, 
LAE077, LAE087, LAE096, LAE101 and LAE102 all 
have relatively small photometrically-estimated REWs 
and upper-limits to their predicted detection significance 
(Section |3] and Table |3j between SL ~ 1 - 6. Further- 
mor e, their fluxes in th e Ly A image are relatively diffuse 
(see lNestor et al.ll2011l Figure 4). We remove these six 
systems from our statistical sample, but are not able to 
rule them out as z ~ 3.09 galaxies with escaping LyC 
flux. 

4.2.4. Summary: LAEs 

Our original s ample of LAE candidates from 
iNestor et all (|201lD contained 27 sources with NB3640 
detections. Our new dataset includes spectra of 26 of 
these 27 sources. We were able to identify an emission 
line that we attribute to Lya in 21 of the 26, although one 
source (LAE034) is at a redshift that is too low for inclu- 
sion in our sample. The other five sources with NB3640 
detections may also be LAEs at z ~ 3.1 as our data 
are not of sufficient sensitivity to detect the expected 
Lya emission line. However, we conservatively remove 
these five systems from our statistical sample. We also 
remove three LAEs with NB3640 detections that show 
evidence for foreground interlopers in the spectroscopic 
data. Thus, our current statistical sample of 91 spectro- 
scopically confirmed z > 3.055 LAEs contains 17 objects 
with NB3640 detections and no evidence for contamina- 
tion of their NB3640 flux by foreground interlopers. We 
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Fig. 17. — The radial surface density of NB3640 detections 
around galaxies in our LBG and LAE samples. The solid red 
histograms include all detected sources; the subset of these 
sources associated with obvious neighbors (which we have ex- 
cluded as possible LyC detections, see i !5.1|) are represented 
by the hatched region. The dashed lines indicate the global 
surface density of sources in our NB3640 magnitude range and 
thus represent the expected levels of contamination, while the 
dotted lines represent the expected la scatter in the contam- 
ination. The excess surface density at low-offsets indicates 
that some of our low-offset LBG NB3640 detections and many 
of the low-offset LAE NB3640 detections are physically asso- 
ciated with the z ~ 3.09 sources and not random foreground 
interlopers. The top panels use the displacement of the LyC 
centroid from that of the R band detection (or BV or NB4980 
for LAEs undetected in R). The bottom panel uses displace- 
ments from the LyA detections of LAEs. As the NB3640 
detections tend to be more spatially coincident with the non- 
ionizing UV continuum than with the LyA emission, the sig- 
nificance of the excess of surface density is greater when using 
the R band offsets. 

stacked the NB3640 images of the 71 LAEs having no 
NB3640 detections and detected no flux down to an 3<r 
limit on the average magnitude of TONB3640 — 28.39. 

5. ACCOUNTING FOR LOWER- REDSHIFT INTERLOPERS 
AND IGM ABSORPTION 



As in INestor et al.l <|2011f ). we use the retained NB3640 
detections in our samples to estimate the contributions 
to the cosmic ionizing background from LBGs and LAEs 
at z ~ 3. However, we must first account for two ef- 
fects. Although we have removed from our samples tar- 
gets with evidence of foreground galaxies contaminat- 
ing their NB3640 flux measurements, the non-detection 
of emission or absorption features belonging to a lower- 
redshift system is not sufficient to rule out the presence 
of an interloper. Such an interloper may, for example, 
lack features strong enough to be detected in our data, 
lie at a redshift such that no strong spectral features fall 
within our wavelength coverage, or not be covered by the 
position of the slit. Not accounting for such interlopers 
would result in an overestimation of cl v c- At the same 
time, the ionizing flux from LyC-leaking galaxies will ex- 
perience an unknown amount of absorption by neutral 
gas in the intervening IGM, which will decrease the ob- 
served value of £LyC relative to the intrinsic value. We 
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account for these two effects in a statistical manner us- 
ing a pair of Monte Carlo simulations, which make use 
of the global surface density of NB3640 detections in the 
relevant magnitude range, and the observed Lya forest 
statistics at z ~ 3. Th e proc edure s are similar to thos e 
used by iShaplev et all (|2006t ) and iNestor et all (|20 1 If ) . 
but employ updated methods and statistics. We briefly 
summarize each method here. 

5.1. Contamination Simulation 

The centroids of the NB3640 detections in our samples 
are generally offset from the corresponding centroids of 
the non-ionizing UV emission. Many of these NB3640 
detections, despite their relatively close proximity on the 
sky to LBGs and LAEs, have clear associations with un- 
related neighboring sources and th us have already bee n 
rejected as possible LyC emission bv lNestor et al.l ((20111) . 
We have now also removed from our samples those galax- 
ies with unconfirmed redshifts, and have rejected NB3640 
detections displaying evidence for contamination in their 
LRIS or NIRSPEC spectra. In Figure [TF] we show the 
resulting surface density of the total (open histogram) 
and rejected (hatched region) number of NB3640 de- 
tections, as a function of offset, for both the LBG and 
LAE samples. For the LAEs, we present surface densi- 
ties computed using the offsets of the NB3640 detections 
from both the UV continuum (i.e., i?-band, or BV for 
LAE081 which is undetected in R) and LyA-hand cen- 
troids. The dashed line in Figure [T7] indicates the global 
surface density of sources in the range of NB3640 mag- 
nitudes spanned by our detections, ps — 0.024 arcsec -2 
over 25 < tonb3640 < 27.25, which corresponds to the ex- 
pected level of contamination. Both the LBG and LAE 
samples have, at relatively small (< 1") offsets, an excess 
of NB3640 detection surface density above the expected 
foreground level. In particular, the excess for the LAEs 
using i?-band offsets is strikingly significant. In order to 
quantify the number NB3640 detections that may be due 
to foreground interlopers, in addition to those already re- 
jected, we performed a Monte Carlo simulation using the 
observed surface densities of sources and predicted levels 
of contamination. 

To account for the increased probability that NB3640 
detections at larger offsets are contaminants, we consid- 
ered each individual detection in turn, computing the 
contamination probability based on offset in the follow- 
ing manner. For the k th detection, we constructed an an- 
nulus, having area A k , encompassing the detection and 
centered on the non-ionizing II V flux centroid. We ap- 
plied this same annulus to each of the N galaxies in the 
sample (N = 41 or 91 for the LBG or LAE samples, 
respectively) and determined the number of previously 
rejected detections, n£ ej -, and putative LyC detections, 
n LyC> contained within the N annuli. We then con- 
structed the probability distribution for the number of 
expected random foreground interlopers, nj ore , in the 
regions spanned by the annuli: 

P(n k fore ) = ( n * re ) p«/- (1 -pjC"-"/".), (i) 

where p = A k x p$ is the probability that a given LBG 
or LAE has a random foreground contaminant in A k 



(see, e.g.. IVanzella et all 120101: INestor et al.l 120111) . In 

each realization of the Monte Carlo simulation, we ran- 
domly chose a value for nh ore from P{n k ore ). Of these 

n fore interlopers, n k e j had already been accounted for. 
Thus the number of additional interlopers predicted to 
lie within the N annuli is n k rrj - n k n . The NB3640 

jot a t-j 

detection in question was flagged as an interloper if 
n fore~ n rej ^ n LyC> an d was retained if n k fore -n k ej < 0. 
Otherwise, we randomly determined if the NB3640 de- 
tection was to be flagged as an interloper based on a 
probability = {n k ore — n k e ^)/n k LyC . We then proceeded 

to the next (i.e., fc th + l) detection and repeated the entire 
process. Once each possible LyC detection in the sam- 
ple had been considered, we recorded the total predicted 
number of additional interlopers. The simulation was 
repeated for a total of 1000 iterations to determine the 
expected number and uncertainty in the average number 
of uncontaminated NB3640 detections. 

The distribution of the number of NB3640 detections 
flagged as interlopers in each realization of the simula- 
tion is shown in Figure [18] for the LBG and LAE sam- 
ples. The widths of the distributions depend slightly on 
the size of annular apertures used in the simulations due 
to NB3640 detections stochastically entering and exiting 
the apertures when their sizes were varied. The vari- 
ability in the distribution widths was small and did not 
depend systematically on the aperture sizes. We used a 
circle of radius equal to the seeing FWHM in the NB3640 
image (0'.'8) for detections with offsets <0'.'4. For detec- 
tions at larger offsets, the radii of the annuli were set such 
that all values of A k were equal. Our simulation suggests 
that 2.6 ± 1.2 of the 6 possible LBGs with LyC detec- 
tions are contaminated by foreground sources. Together 
with the three LBGs showing evidence for contamination 
in their spectra discussed in Section 14.1.21 the resulting 
contamination rate is 62 ± 13% (i.e., 5.6±1.2 of 9). The 
contamination-corrected detection rate for the sample as 
a whole is 8 ± 3% (3.4 ± 1.2 of 41). 

Using i?-band offsets, our simulation suggests that 
3.8± 1.3 of the 17 possible LAEs with LyC detections are 
contaminated. Together with the three LAEs with con- 
tamination discussed in Section 14.2.21 the resulting con- 
tamination rate for the LAE sample is 34± 7% (6.8± 1.3 
of 20) and the contamination-corrected detection rate is 
15 ± 1% (13.2 ± 1.3 of 91). If we instead use LyA off- 
sets, the predicted number additional contaminants be- 
comes 6.4±1.9, the contamination rate becomes 47±10% 
(9.4 ± 1.9 of 20) and the detection rate becomes 12 ± 2% 
(10.6 ± 1.9 of 91). 

In addition to estimating the number of foreground in- 
terlopers, our simulation computes the (contamination- 
corrected) average LyC and R magnitudes and uncer- 
tainties. The uncertainties include sample variance com- 
puted by first randomly reassigning individual magni- 
tudes based on the measured magnitude and error, as- 
suming Gaussian magnitude uncertainties determined 
from our photometric simulations, and then bootstrap 
re-sampling each data set. When computing the sample- 
average LyC magnitudes, we alternately assumed that 
NB4630 non-detections and detections flagged as inter- 
lopers had flux levels equal to zero or the maximum aver- 
age magnitude consistent with the (lcr) limits set by our 
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Fig. 18. — The distribution of the number of contaminated 
NB3640 detections as predicted by our Monte Carlo simu- 
lations, excluding the three LBGs and three LAEs showing 
evidence for foreground contamination in their spectra. The 
upper panel shows the results for the LBGs, which have a to- 
tal of six possible LyC detections, while the lower panel shows 
the results for the LAEs, which have a total of 17 possible 
LyC detections. In the simulations, NB3640 detections with 
larger spatial offsets from the corresponding non-ionizing UV 
flux centroids are more likely to be flagged as interlopers. For 
the LAE sample, we considered offsets from both the R band 
(black solid histogram) and LyA (red dashed histogram) flux 
centroids. 

stacking analysis. The resulting contaminated-corrected 
sample-average NB3640— R colors are listed in Table SJ 

5.2. IGM Simulation 

To investigate the attenuation of escaping LyC flux by 
absorption from the IGM, we used the observed column 
density and frequency distributions of the Lya forest over 
the relevant redshift range, 1.7 < z < z source , to model 
500 random sightlines through the IGM for each LBG 
and LAE redshift. We then computed the attenuation 
of the continuum flux in the NB3640 filter due to the 
randomly generated neutral clouds. In this manner we 
determined the fraction of flux transmitted, P , for each 
of the 500 si ghtlines. The proce ss was identical to that 
described in iNestor et al.l (|201lD except that we made 
use of updated estimates of the IGM opacity (Rudie et 



al. 2012, in prep.) and a more precise treatment of the 
higher-order Lyman absorption lines from each cloud. As 
the galaxies in our sample lie either in or behind the 
SSA22a protocluster, it may be expected that the IGM 
in the proximate foreground of our sources is either more 
opaque due to the mass-overdensity of the region, or less 
opaque due to the overdensity of ionizing emissivity. In 
either scenario, the difference in average IGM opacity 
should manifest in differences in rest-frame UV colors 
that span Lya emission and/or the Lyman limit. We 
compared the U — G and G — R colors for our sources 
at 3.06 < z < 3.12 with those of field LBGs in the same 
redshift range and found no statistical differences in their 
distributions, suggesting that any such effect is small in 
relation to our other uncertainties. 

For each of the N redshifts in a sample, the expec- 
tation value for the fraction of transmitted flux, (t(z)), 
is simply the average of the 500 simulated t°{z) values. 
For a single galaxy at z ~ 3.09, the distribution of t 
is very broad, ranging from w — 60% (see, e.g., Fig- 
ure 8 of INestor et al.l 1201 ID . The uncertainty in t(z) is 
thus a major uncertainty in our estimate of the contri- 
bution to CLyC from each individual source. The LBG 
and LAE sample- average LyC flux values, however, can 
be corrected with much less uncertainty. We define the 
sample-average transmission as 

^sample — N i 

where Fi are the N individual intrinsic (i.e., prior to 
any IGM absorption) NB3640 fluxes, and t(zi) are the 
actual IGM transmission values. As before, N = 41 or 
91 for the LBG or LAE samples, respectively. Although 
both the individual Fi and t{zi) values are unknown, as 
the two sets are independent the expectation value for 
isample is simply the average of the N values of (t(z)}, 
and is thus independent of the values of Fi and t(zj). The 
uncertainty in i sample, o-j, does depend on the t(z) and F 
distributions, however. While we have an accurate model 
for the probability distributions for the t(zj) values, th e 
Fi values are poorly constrained. In INestor et al.l (|2011l) . 
we estimated ai by randomly drawing the N values of 
t(zi) from the corresponding t 3 (zi) values to compute 
isample- We computed i sa mpie in this manner 1000 times, 
and equated o~i to the standard deviation of the resulting 
isample values. That procedure is equivalent to assuming 
the ./V values of Fi are all equal, and will underestimate 
o~i for more realistic distributions of F. 

To improve upon our estimation of o~i , we assumed an 
exponentially decreasing function for the intrinsic flux 
probability distributions: p(F) oc e~ F ^ . This choice of 
parametrization has the advantage that the results are 
only mildly sensitive to the e-folding parameter, (3. To 
determine the best choice for /3, we convolved p{F) with 
the t(z) distributions for a range of (3 values. We then 
used the resulting attenuated-flux probability distribu- 
tions to compute the likelihood of our contamination- 
corrected LBG and LAE NB3640 data sets, retaining the 
value of j3 that maximized these likelihoods. To deter- 
mine o~i j we randomly selected the N Fi values from the 
maximum likelihood exponential distributions. For each 
flux value, we also randomly selected t(zi) from one of 
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the sets of P(zi) values. These F{ and t(zi) values were 
used with Equation [2] to determine t sam pie- This process 
was repeated 1000 times, and ui was set equal to the 
standard deviation in the 1000 simulated i sam pi e values. 

We found that an exponential form for p(F), when 
convolved with our modeled IGM transmission values, 
resulted in a qualitatively acceptable match to our ob- 
served NB3640 fluxes. However, as the actual F distri- 
butions are only poorly constrained by the data, it is 
likely that other functional forms are also able to match 
the observations and may result in different estimates of 
<7f. In particular, if the true p(F) distributions contain a 
spike at zero flux (as suggested by our non-detection of 
NB360 flux in the stacked images of the non-detections), 
we would still be underestimating aj. Nonetheless, as 
our exponential model represents an improvement over 
the past method which effectively assumed p(F) = con- 
stant, and the error budget in thyC is not dominated by 
the uncertainty in t samp i e , we adopt the o~i values deter- 
mined with an exponential p(F). 

For our LBG sample, we found tsampie = 0.298 ±0.040, 
and for our LAE sample we found t samp i e = 0.320±0.027. 
We used these values to correct the sample-average 
NB3640 magnitudes and NB3640-.R colors. These col- 
ors can be expressed in terms of the sample-average 
UV-to-LyC flux-density ratios, rj = {F uv /F Ly c}- Af- 
ter applying the contamination and IGM corrections, 



we find t\lbg = 18.0 



+34.8 
-7.4 



for our LBG sample, and 



Vlae — for our LAE sample. These values are 

consistent with, though larger than, those estimated by 
INestor et all (IMl : T] LB g = ll^t"^ 3 for LBGs and 
Vlae = 2.2+p'g for LAEs. The relative uncertainties 
in our updated estimates of 77 are larger than those in 
our previous estimates, due to an improved treatment 
of the errors. The updated samples also contribute to 
the increased relative uncertainties, as the newly con- 
firmed LBGs have, on average, larger i?-band photomet- 
ric uncertainties than those in the previous sample, and 
the spectroscopically-confirmed LAE sample used here is 
smaller than the full photometrically-selected LAE sam- 
ple. Nonetheless, by refining our samples and including 
empirical evidence for the presence or absence of fore- 
ground contamination in individual systems, our revised 
values represent improved estimates of 77 for LBGs and 
LAEs at z ~ 3. We list the colors and r] values in Tabled 
for our raw LBG and LAE samples, for those samples 
after application of the contamination corrections, and 
after correcting for both contamination and IGM absorp- 
tion. The contamination- and IGM-corrected values are 
used below to estimate tLyC- 

6. RESULTS 

6.1. Revised Estimate of z ~ 3.01 LyC Emissivity 

One of the primary goals of this paper is to deter- 
mine the global luminosity space density of ionizing ra- 
diation, ehyCi contributed by star- forming galaxies at 
z ~ 3. As the SSA22 protocluster represents an over- 
dense volume of the Universe, computing tLyC directly 
from our data would overestimate the space density of 
ionizing flux. However, the UV-continuum luminos- 
ity functions of LBGs and LAEs have been measured 
over relatively representative volumes, at rest-frame ef- 
fective wavelengths close to those of our i?-band data 
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Fig. 19. — The rest-frame A ~ 1600 A absolute magnitude 
distribution for our LBG (solid black) and LAE (dash-dot 
red) samples. The hatched regions indicate galaxies with 
NB3640 detections. The faintest bin in the LAE histogram 
also contains the 10 LAEs fainter than our detection limit of 
Mab = —18.3. The top axis indicates the corresponding lu- 
minosity relative to the LBG characteristic luminosity L* LBG . 
Our LBG sample is dominated by galaxies brighter than 
Mab = —20, corresponding to L ~ 0A6L LBG , while our LAE 
sample is dominated by sources fainter than Mab = —20. 



A ~ 1600A), bv iReddv etaTl (pool and lOuchi et~aH 
2008), respectively. Thus, we can use the values de- 
termined for the sample-average UV-to-LyC flux-density 
ratios together with the corresponding luminosity func- 
tions to determine £LyC : 



<L . :l r I ^L$dL, 



(3) 



where L refers to the nora-ionizing UV continuum lumi- 
nosity, and r\ may be a function of L. When determined 
in this fashion, the value computed for thyC will depend 
on the choices of shape and normalization for the lumi- 
nosity function, and the range in luminosity over which 
Equation [3] is in tegrat ed. 

IReddv et all (l200l give Schechter function parame- 
ters for the A ~ 1700 A UV-continuum luminosity func- 
tion of LBGs at z ~ 3 of 4>* LBG = 1.66 x 10~ 3 ™— ~ 3 



M 



Mpc" 



LBG 



-20.84, and a LB G = -1.57. For LAEs at 



z ~ 3.1. lOuchi et"aH l|20Q8h determine A - 1600A UV- 
continuum luminosity function parameters of (\>\ AB — 



0.56 x Mpc~ 



M 



LAE 



-19. 8, and a^AE 



— 1.6. However, their luminosity function is determined 
for LAEs having REW > 64A, while our sample has 
R EW > 20 A. From the REW distribution determined 
bv lGronwall et all (|2007f ), we estimate that LAEs repre- 
sented by our sample are a facto r of » 1.8 more com mon 
at z ~ 3 than those described bv lOuchi et all (|2008[ ). We 
therefore adopt (j)* LAE = 1.01 X 10" 3 Mpc -3 for our LAE 
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sample 

The absolute magnitude distributions for our LBG and 
LAE samples are shown in Figure [T9l Spectroscopic con- 
firmation of color-selected LBGs in the SSA22a field was 
limited to 72 < 25.5. Consequently, the bulk of our LBG 
sample is brighter than Mab = — 20. In contrast, our 
LAE sample, which includes galaxies down to our photo- 
metric limit of R = 27.3 (or M AB ~ -18.3 at z ~ 3.09), 
is dominated by galaxies fainter than Mab — —20. Thus, 
we first compute the contributions to CLyC from LBGs 
and LAEs separately, and only from sources within the 
UV-continuum luminosity ranges over which we deter- 
mined the values of ry: Mab < —20.0 or our LBG sam- 
ple, and —20 < Mab < —18.3 for our LAE sample. 
The absolute magnitude limits Mab = —20.0 and —18.3 
correspond to L m i n — QAQL* LBG and Q.\L* LBG , respec- 
tively. With these magnitude ranges, the parameters 
for the respective luminosity functions discussed above, 
and riLBG = 18.0t? 4 4 8 for LBGs and j] LA e = 3.7±^ 
for LAEs as determined in Section [5J we find £ L yc = 
5.2±|;2 x 10 24 ergs s" 1 Hz" 1 Mpc~ 3 for M AB < -20.0 



and e LAE 



6.7+3:4 x 10 24 ergs s" 1 Hz" 1 Mpc" 3 for 



+3.6 

LyC ■ 

-20 < M AB < -18.3. 

We now turn our attention to the total contribution 
to CLyc from all star-forming galaxies brighter than 
Mab — —18.3 at z ~ 3. First, however, we must address 
the principle source of the large (factor of ~ 5) differ- 
ence between the values determined for rj in our LBG 
and LAE samples. As the LBG and LAE samples are 
dominated by galaxies with relatively bright and faint 
UV-continuum luminosities, respectively, it may be that 
rj has a strong dependence on galaxy luminosity. In this 
scenario, which we refer to as our luminosity-dependent 
r\ model, we identify t\lbg with all star forming galaxies 
having Mab < —20.0 and and vlae with all star forming 
galaxies having —20 < Mab < —18.3. Here w e make the 
approximation that the iReddv et al.l (|2008[ ) luminosity 
function, which is determined from color-selected galax- 
ies, is a complete census of all star-forming galaxies with 
Mab < —18.3; i.e., that there are no sources brighter 
than this limit that contribute to £LyC at z ~ 3 with UV 
colors such that they would not be selected as LBGs. 
With these assumptions we can integrate Equation [3] 
piecewise: 



lurn. — dep. 

C ~LyC 



1 



0.46L* 



VLAE Jq.1L' 
1 



Vlbg J0A6L 



L§ LBG dL 



L<Z> LBG dL 



(4) 



finding e ^- de f- = 32.2™ ergs s" 1 Hz" 1 Mpc" 3 , 
with contributions from sources fainter and brighter than 
M AB = -20 of 27.0± 44 : 4 x 10 24 ergs s" 1 Hz" 1 Mpc~ 3 
and 5.2l 3 4 x 10 24 ergs s -1 Hz -1 Mpc -3 , respectively. 
This value of e^~ dep - is » 60% of that estimated in 

4 Although the UV luminosity function of LAEs having 20A < 
REW < 64A (which compose ~ half of LAEs with REW > 20A) 
may, in principle, have different values of M* and a compared to 
LAEs with REW > 64A, we proceed under the assumption that 
the Ouchi ct al. (2008) values for M* and a apply for our entire 
range of REW. 



iNestor et al.l (|201 lh . primarily due to our refined values 
for t)lbg and t)lae- 

Alternatively, the large difference in rj between the 
LBG and LAE samples may be driven by some prop- 
erty of star-forming galaxies associated with LAEs. In 
this scenario, which we refer to as our LAE-dependent 
r\ model, we can estimate the total value of thyc by 
assuming r\LAE holds for all LAEs over the full range 
Mab < —18.3, and t]lbg is representative of all LBGs 
with Mab < —18.3 that are not also LAEs. The frac- 
tion of LBGs that are not LAEs is determined from the 
LBG and LAE space densities, which can be computed 
by integrating the respective LBG and LAE luminosity 
functions, p = . $eLL0 Setting L min = 0.1L LBG , 

we estimate that 23% of LBGs are also LAEs with REW 
> 20A, which is consistent with past results at z ~ 
3 (iSteidel et al.l [20001 IShaplev et al][200l iKornei et al.l 
2010). Here we have made the approximation that LAEs 
comprise a sub-sample of LBGs; i.e., that all LAEs at 
z ~ 3 with Mab < —18.3 wou l d mee t the color-selection 
criteria used by IReddv et al.l ((2008) . Equation [3] then 
becomes: 



LAE- 
€ LyC 



0.77 

Vlbg J i l 
1 



LBG 



VLAE Jo.1L 



LAE 



dL 



dL, 



(5) 



LAE— dep. 
e LyC 

Hz- 1 



16.8±g-J 



resulting ir 

10 24 ergs s _1 Hz -1 Mpc -3 , with contribu- 
tions from non-LAE and LAE galaxies of 



,+5.1 



-5.5 



10 2 



ergs 



Hz~ 



8.6i^ x 10 24 ergs s" 1 Hz" 1 Mpc -3 



Mpc -3 and 
respectively. 

This value of ^yc~ dep ' is similar to that estimated 
for just LBGs with L > 0.1L* in INestor et"aH (|2011h . 
and ~ 30% of the total value of tLyC estimated by 
INestor et all (|201lD . This difference is due to the 
combination of our improved measurements of t]lbg 
and tjlae and the reduced contribution from sources 
with low values of 77 in our LAE-dependent model. 

It is difficult, given the present data, to distinguish 
which of the luminosity-dependen t or LAE-dependen t 
rj models is more appropriate. In N estor et al.l (201 1). 
we investigated differences in the average Lya emission 
strengths of systems with NB3640 detections compared 
to those without detections by creating stacks of LRIS 
spectra of our LBGs, and by comparing the distributions 
of photometrically-estimated Lya REWs of our LAEs. 
For both LBGs and LAEs, we found that systems with 
NB3640 detections tend to have weaker Lya emission 
relative to those not detected in NB3640. We repeated 
these tests using our refined LBG and LAE samples and 
excluding systems showing evidence for contamination 
in their spectra, and found trends consistent with our 
previous results. Thus, it is unlikely that Lya emission 
strength has a direct influence on the average value of 77. 
However, the LyC properties of z ~ 3 galaxies may be de- 
pendent on some other property or combination of prop- 
erties, such as star formation surface density, stellar pop- 

In this simple model, we do not account for the observed lu- 
minosity dependence of the fraction of LBGs that are also LAEs 
(see, e.g. JStark et~aLll20TOh . 
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ulation age, metallicity, etc., that is more typical of the 
LAEs in our sample compared to our LBGs, independent 
of UV-continuum lu minosity. For example, in a sample 
of local LBG analogs. iHeckman et al.l (|2011[ ) present evi- 
dence for small covering factors for optically thick neutral 
gas in the most compact galaxies in their sample. Indeed, 
in the available HST/ ACS F814W imaging, the LAEs in 
our sample do appear, qualitatively, more compact on av- 
erage than the LBGs. Nonetheless, as our LBG and LAE 
samples are largely distinct in magnitude range, future 
work including including fainter LBGs and larger LAE 
samples is needed to clearly differentiate between the two 
proposed scenarios. 

We summarize the various contributions to ^LyC in 
Table [5] The uncertainties in the values of thyC are 
dominated by the uncertainties in our estimates "q, 
which are large compared to the errors in the first mo- 
ments of the luminosity functions. The uncertainties in 
€l v c do not, however, include systematic errors arising 
from our choice of lum inosity functions. For example, 
iReddv fc Steidell ([2009) find a steeper faint-end slope at 
z ~ 3 previous works (see discussion in IReddv fc Steidel 
2009). Furthermore, in estimating r\ we have only ac- 
counted for the contribution from galaxies brighter than 
0.1L*, i.e., the luminosity range over which we have em- 
pirica l estimates of 77. We note that, for the IReddv et al.l 
(2008) luminosity function, galaxies fainter than 0.1L* 
contribute ~ 40% of the luminosity density at 1600A. 

Given the spectral shape of the ionizing continuum flux 
and the mean free path to ionizing radiation in the IGM, 
A m fp, the value of tLyC implies a corresponding (proper) 
hydrogen photoionization rate in the IGM, Thi- If we 
assume a power law form for the LyC flux density such 
that f v oc is~ a , then 



HI 



(l + z) 3 cr H iA m f p 
h (3 + a) 



all, L>0.1L* 
-LyC 



(6) 



where h is Planck's constant, and oh i = 6.3 x 10 -18 cm~ 2 
is the atomic hydrogen photoi onization cros s secti on at 
the Lyman limit. Following INestor et al.l (|2011[ ). we 
adopt a — 3 and A m f p = 75.6 Mpc. For the two es- 



timates of ^LyC determined above. 

2.7iJ' g x 10~ 12 s _1 in our luminosity-dependent 77 model, 



n 1 T^lum. — c 

we find I H j 



and r£f E ~ dep - = 1.4 ± 0.6 x 10~ 12 s" 1 in our LAE- 
dependent 77 model0 We list these values of Thi in Ta- 
bled tog ether with estimat e s of T hi in the z ~ 3.1 Lyq 
forest by iMeiksin fc White! (12004ft [Bolton fc Haehnelti 
(|2007fl . and iFaucher-Giguere et"atl |2008|). The values 
we obtain for Thi are larger than the values reported 
in the literature by a factor of « 2 — 4, although they 
are consistent given the uncertainties with the literature 
values at ~ la and ~ ler, for the luminosity- and LAE- 
dependent 77 assumptions, respectively. 

6.2. Escape Fractions 



6 More recent estimates of the mean free path of ionizing photons 
at z ~ 3 suggest A m f p = 63 Mpc in the IGM, and A m f p = 49 Mpc 
if the opacity contributed by the circum-galactic medium of LBGs 
is included (G. Rudie, private communication). Adopting these 
values would lead to estimates of Thi that are 83% and 65%, re- 
spectively, of the values presented here. 



It is often of interest to estimate the fraction of LyC 
radiation produced by star formation that escapes into 
the IGM, f^yc. This value can be useful for constraining 
models of star formation histories, star formation feed- 
back in the interstellar medium (ISM), etc. Additionally, 
/esc can provide a useful parametrization in reionization 
models. The value of f^Jj, can be estimated empirically 
by 

-stars 



/•LyC = 
J esc 



V 



f, 



uv 

esc ' 



(?) 



where r7 stars is the intrinsic ratio of UV-to-LyC luminos- 
ity densities in star-forming regions, and /HY is the es- 
cape f raction of non- ionizing UV radiation. IReddv et al.l 
(2008) have estimated « 20% in LBGs at z ~ 

3. LAEs, however, exhibit bluer average UV continua 
and smaller UV attenuation, with .ff; 30% (e.g ., 
lOuchi et al.ll200l iKornei et aJl[20lot iBlanc et al.|[20Tlh . 
The term r] stars /r], which is often referred to as the rela- 
tive escape fraction, is a measure of the attenuation of the 
UV flux relative to that of the LyC flux. In this work 
we have estimated r\ for samples of LBGs and LAEs. 
The value rj stars , however, must be obtained from spec- 
tral synthesis models, and depends on the ages of the 
stellar populations. Even for a given star-formation his- 
tory, various models produce differing values for rj stars . 
To illustrate the dependencies of the inferred value of 
feJc C on different models and parameters, we determined 
stars ugm g the st e llar p opulation synthesis models of 
IBruzual fc Charlotl (12003 . BC03), as well as a recent set 
of model integrated spectral energy distributions pro- 
duced by The Binary Population and Sp ectral Synthesis 
code (BPASS . lEldridge fc Stanwayl2009t ) . which includes 
the effects of massive binary stars and nebular emission. 
For each of the models we adopted a constant SFR and 
two different metallicites (Z = 0.004 and 0.020) and mea- 
sured rj stars , using our NB3640 and R filter passbands, at 
three ages (T = 10 6 , 10 7 and 10 8 yr). The model values 
we obtain for rj stars range from ~ 1.3 for the T = 10 6 yr, 
Z = 0.004 population in the BPASS models, to ~ 6.4 
for the T = 10 8 yr, Z = 0.02 population in the BC03 
models. The complete set of modeled rj stars values are 
given in Table [71 

For a given combination of ?7 and f^, the modeled 
values of rj stars imply a corresponding f^ c . We show 
these fe S y c c values in Table [7J for each value of rj stars , 
assuming t\lbg — 18-0 for LBGs and tj^ae = 3.7 for 
LAEs, and adopting, in turn, f%%= 0.2, fg%= 0.3, and 
the limiting dust-free case where 1- The f^ s v c c 

values vary by factors of ~ 5 — 7 for a given rj and 
fesc- They vary from 1% for LBGs in the youngest 
JYZ= 0-2 BPASS model, to - 50% for LAEs in the 
oldest fesY = 0-3 BC03 models, and reach greater than 
unity for several of the dust-free LAE models. For a given 
/^Y an d T), fes V c C varies most with stellar population age, 
although the choice of BC03 or BPASS model also affects 
the inferred value of f^ s y c c by as much as a factor of ~ 2. 
Over the range investigated, metallicity has only a small 
effect on f^ s v c C - There are other model inputs that we 
have not varied, such as the stellar initial mass function 
or star-formation history, that are also likely to affect 
■q stars and therefore the determination of fe S v c C - With 
these caveats in mind, LBGs at z ~ 3 have a median 
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age of ~ 300 Myr dKornei et al.ll2010D wi t h metallicities 
approaching ~ solar (|Pettini et al.l 120011 : iShaplev et al.l 
12004 IMaiolino et"all I2008D and fg^- 0.2, suggesting 
fLyCLBG _ 5 _ ?% _ High-red shift LAEs, in con - 
trast, have typic al ages ~ 20 Myr dGawiser et al][ 2007'). 
low metallicity (jFinkelstein et all 120111 : iNakajima et al.l 
12011 . and f™~ 0-3, suggesting /i^ C ' LAE ~ 10 - 30%. 
In general, we caution against the naive use of LyC es- 
cape fractions in reionization models without considera- 
tion of the values of rj stars and used in their deter- 
mination. 

In Table @] we also list the values of r] for only galax- 
ies with LyC detections. Our sample of LBGs with LyC 
detections has r/ = 1.7l 'g « ^stars mos ^ f our moc i_ 
eled j] stars values, with the exception of the older BC03 
models, implying a ~ unity relative escape fraction in 
LBGs with LyC detections. Note that, as LyC radia- 
tion is more susceptible than non-ionizing UV radiation 
to attenuation by gas and dust, rj ss jj stars implies little 
attenuation and therefore a large value of /J^f. For the 
LAEs with NB3640 detections, our IGM-corrected value 
T) = 0.41° i is inconsistent with all of our modeled val- 
ues of r) stars at more than 3<r (note that it is unphysical 
to have r\ > rj stars ). We can bring our LAE detections- 
only value for r\ in line with the model predictions if we 
assume that we can only detect LyC in our LAE sam- 
ple along fortuitously clear IGM sightlines, and therefore 
can neglect the IGM correction, resulting in -q = 1.4^4. 
This scenario, however, is very unlikelyQ Deciphering 
the correct interpretation of these surprising low values 
of rj is a primary goal of our ongoing work, which in- 
cludes H ST /UVTS LyC and UV imaging of many of our 
NB3640 detections. Better multi-wavelength constraints 
on the stellar populations of NB3640-detected LAEs will 
also be a key component of determining their underlying 
nature. 

In Section [Ol we found a LyC detection rate of ~ 8% 
in our LBG sample, and ~ 12 — 15% in our LAE sample. 
While it is not clear if the LyC properties of a "typi- 
cal" z ~ 3 LBG or LAE are similar to that of the aver- 
age system, these rates do indicate the solid angle over 
which LyC radiation escapes, at a level above our de- 
te ction limit, average d over all galaxies in each sample. 
In lNestor et al.l (|2011| ). we proposed a "blow out" model 
in which feedback from regions of dense star formation 
clear portions of the ISM of gas and dust. When viewed 
along favorable sightlines, such regions appear to have 
large escape fractions, as we find in our LBGs and LAEs 
with LyC detections. Galaxies that either have failed to 
sufficiently remove their ISM over significant solid angle, 
or are viewed along unfavorable sightlines, will appear 
to have negligible escape fractions. The strong LyC- 
flux upper-limits in systems without NB3640 detections, 
derived from our stacking analysis, are consistent with 
this picture. Additionally, in the subsample of our LAEs 
having HST imaging[3 we find no significant difference 

7 The average of the largest 12% (equal to the LAE LyC de- 
tection rate) of our modeled IGM transmissions is 0.566. Even 
adopting this value for the IGM correction we find r\ = 0.79, which 
is still significantly below our lowest predicted value of •q stars . 

8 The number of LBGs with NB3640 detections and HST imag- 
ing is too small to make meaningful comparisons to the sample of 



in the distributions of sizes or surface brightnesses for 
sources with and without NB3640 detections. The simi- 
larity of these properties between NB3640 detected and 
non-detected LAEs is consistent with a scenario in which 
viewing angle is a significant factor in the ability to de- 
tect escaping LyC. 

7. DISCUSSION AND SUMMARY 

Directly studying galaxies in the ionizing continuum is 
a difficult endeavor. The dearth of QSOs above z ~ 3 
suggests that, at high redshift, the ionization balance in 
the IGM is maintained by LyC flux escaping from star- 
forming galaxies. However, the increasingly opaque IGM 
makes the detection of a ny escaping LyC flux unlikely 
above z ~ 3.5 (see, e.g., IVanzella et al.ll2012f) . Below 
z ~ 2.4, the redshifted Lyman limit falls below the atmo- 
spheric cutoff requiring observations from space to detect 
ionizing flux. Current observations at z m 1.3 sampling 
rest- fra me A ~ 700 A have result ed only in upper limits 
to e LyC (jSiana et alj I2007L 12010). The non-detection of 
LyC emission at z ~ 1 together with the apparent need 
for a galaxy contribution to €LyC at high redshift implies 
that clv C evolves strongly over z ~ 1 — 3 (se e, e.g., Fig- 
ure 9 of lNestor et"aLll20lU Unoue et al.H2006f) . We note, 
however, that if the average UV to LyC flux-density ratio 
rj is luminosity-dependent, the non-detection of LyC in 
star forming galaxies at z ~ 1.3 could be d ue in part to 
a selec tion bias, as the galaxies observed by iSiana et al.1 
(|2010t) are preferentially bright. 

In any case, short of a large space-based program, the 
window 2.5 < z < 3.5 offers the best opportunity to 
identify and study LyC-leaking galaxies. As this and 
other recent works have shown, however, difficulties re- 
main even at z ~ 3. While the IGM is more transparent 
to ionizing radiation at z ~ 3 compared to higher red- 
shift, there remains an unknown level of IGM attenuation 
of escaping LyC flux. The only practical way to correct 
for the attenuation is by applying a statistical correction, 
determined through modeling of the IGM, to the sample- 
average LyC properties. In this work, we have improved 
upon past efforts to estimate this correction by using an 
updated estimate of the IGM opacity at z ~ 3, includ- 
ing a more precise treatment of the higher-order Lyman 
absorption lines, and developing a new method to obtain 
a more realistic estimate of the transmission uncertain- 
ties. For the sample-average IGM correction method to 
be effective, we require a relatively large statistical sam- 
ple. This necessity presents a second difficulty, as each 
galaxy requires a deep observation at a large observa- 
tory to be detected in the UV and LyC at these red- 
shifts. Thus, we conducted our observations in a single 
field that contains a large overdensity of galaxies in a 
narrow range of redshift, using Keck and archival Sub- 
aru and HST imaging data, and new Keck/LRIS multi- 
object spectroscopy. Our resulting LBG and LAE sam- 
ples are the largest of their types to have been imaged 
below the Lyman limit in an attempt to identify LyC 
flux. T he samples presente d here are improvements of 
those in iNestor et"atl (|2011l ). through a ~ 60% increase 
in the number of LBGs and the spectroscopic confirma- 
tion of > 90% of our LAE candidates. Finally, perhaps 
the most critical impediment to unambiguous identifi- 

LBGs without NB3640 detections. 
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cation of leaking LyC at z ~ 3 is the high sky surface 
density of foreground sources with observed-frame blue 
magnitudes similar to expected z ~ 3 LyC magnitudes. 
Our new spectra mitigate this problem by allowing us 
to directly search for spectral features from galaxies at 
lower redshift nearby on the sky (and thus fal l ing on the 
slit) to the z ~ 3 galaxies. As in lNestor et al.l (|201lD . we 
also statistically correct the samples for yet unidentified 
foreground contaminants, here using a slightly refined 
method. Overall, in addition to clarifying the interpre- 
tation of the NB3640 detections of several specific LBGs 
and LAEs, our new data have led to improved estimates 
of the galaxy contributions to the ionizing background 
flux and ionization rates. Our main results include: 

1. We obtained spectra of 41 LBG candidates and 
confirm redshifts z > 2.45 for 26. Of these 26, 
16 have z > 3.055 such that our NB3640 filter 
is opaque above the redshifted Lyman break. To- 
gether with the 25 previously identified z > 3.055 
LBGs in the field, our NB3640 image samples be- 
low the Lyman limit for a total of 41 confirmed 
LBGs at z > 3.055. 

2. We obtaine d spectra of 9 6 of t he 110 LAE candi- 
dates from I Nest or et al.l (| 2 1 If) . We confirm red- 
shifts 3.057 < z < 3.108 for 87 of the candidates. 
An additional four candidates, which are also part 
of our LBG sample, have previously determined 
redshifts z > 3.055. Thus, our our NB3640 im- 
age samples below the Lyman limit for a total of 
91 confirmed LAEs. 

3. Nine galaxies in our LBG sample are detected in 
NB3640. The spectra of three of these detections 
contain evidence for the presence of a lower-redshift 
interloper, and thus are removed from consider- 
ation as LyC leakers. Our Monte Carlo simula- 
tion suggests that an additional 2.6 ± 1.2 of our 
NB3640 detections are contaminated. The result- 
ing predicted contamination rate is 62 ± 13%, and 
the predicted LyC detection rate is 8 ± 3%. 

4. We detect 20 galaxies in our LAE sample in the 
NB3640 image. The spectra of three of these detec- 
tions contain evidence for the presence of a lower- 
redshift interloper, and thus are removed from con- 
sideration as LyC leakers. Our Monte Carlo simu- 
lation suggests that an additional 6.4 ± 1.9 of our 
NB3640 detections are contaminated. The result- 
ing predicted contamination rate is 47 ± 10%, and 
the predicted LyC detection rate is 12 ± 2%. 

5. Using the sample- average NB3640— R colors, we 
determined the sample-average UV-to-LyC flux- 
density ratios r\ — 18.0^y 4 4 8 for our LBG sample 

and r\ = 2>.T^\ for our LAE sample. We then 
used these ratios, together with the z ~ 3 LBG 
and LAE luminosity functions to estimate the con- 
tributions to thyC from each sample, over the lumi- 
nosity ranges for which r\ was constrained: t^yS = 
5.2j^xl0 24 ergss- 1 Hz" 1 Mpc~ 3 for Mab < -20 
and e£$f = 6.7±|;^ x 10 24 ergs s" 1 Hz" 1 Mpc" 3 
for -20 < Mab < -18.3. The total 



galaxy contributions are e l ^^~ de P- — 32.2^} 2 '4 x 
10 24 ergs s _1 Hz -1 Mpc -3 in our luminosity- 
dependent r] model or e I j^yQ~ dep ' — 16.8+g'5 x 
10 24 ergs s _1 Hz -1 Mpc -3 in our LAE-dependent 
rj model, for galaxies with Mab < —18.3 which 
corresponds to L ~ Q.\L* LBG at z ~ 3. 

6. The intergalactic hydrogen photoionization rates 
L inferred by our values of CL y c, L^'J™' dep ' — 

x 10 -12 s _1 in our luminosity-dependent r] 

model, and r^ E ~ dep - = 1.4 ± 0.6 x 10~ 12 s" 1 in 
our LAE-dependent 7/ model, are larger than but 
consistent with published values for L measured in 
the Lya forest. 

7. We use two suites of stellar population synthesis 
models to predict the intrinsic UV-to-LyC flux- 
density ratio in star-forming regions, ■q stars , find- 
ing values that range from 1.3 — 6.4. The inferred 
sample-average values of for our samples de- 
pend on the assumed values of rf tars and jYsc ■ 
Despite the uncertainty in these two values, our 
best estimates of f^ s v c c are ~ 5 — 7% for LBGs and 
- 10 - 30% for LAEs. In our "blow out" model 
for LyC escape, we predict near-unity apparent es- 
cape fractions in galaxies with LyC detections, and 
negligible escape fractions in other systems. 

Despite our significant progress, several outstanding 
issues remain. In order to assess the robustness of our 
contamination corrections, it is important to determine 
unambiguously the LyC or interloper nature of as many 
of our NB3640 detections as possible. For many of our 
galaxies with NB3640 detections, our LRIS and/or NIR- 
SPEC data enable us to associate a redshift with one or 
more specific regions (or "clumps" ) of non-ionizing UV 
flux. With the data at hand, we are already able to place 
high confidence on the association of NB3640 detections 
with z > 3.055 clumps in at least two of our objects, 
MD46 and LAE053. We expect to increase the number 
of such high-confidence LyC identifications in the near 
future with our ongoing programs. Notably, the resolu- 
tion afforded by our ongoing H ST /WIS cycle 19 imag- 
ing program will allow the measurement of Fjjv I FhyC m 
the individual clumps, thus providing strong constraints 
on the nature of the NB3640 detections, i^-band and/or 
Spitzer /IRAC data, together with our rest-frame UV ob- 
servations, will enable investigations of the stellar popu- 
lations of galaxies with and without LyC detections. It is 
also necessary to extend our investigation to other fields 
in order to ensure that our results are not a phenomenon 
(or, e.g., a systematic bias) unique to the SSA22a field. 
To this end, we are in the process of finalizing a similar 
study in an independent field at z ~ 2.85 (Mostardi et 
al., in prep). Finally, the recently commissioned near-IR 
multi-object spectrograph MOSFIRE at the Keck obser- 
vatory will allow the study of rest-frame optical spectral 
features of galaxies with LyC detections. 
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APPENDIX 
LAE LYMAN-CONTINUUM CANDIDATES 

LAE010 

We show the H ST /ACS imaging and one-dimensional shallow-mask spectrum of LAE010 in the top panel of Figure[7l 
An emission line is clearly detected at A = 4979A, indicating a redshift of z = 3.096. The shallow-mask slit also covers 
the position of the NB3640 detection. No additional emission features are detected in the spectrum. Thus, we retain 
LAE010 as a possible LyC-leaking galaxy 

LAE016 

We show the HST/ ACS imaging and one-dimensional shallow-mask spectrum of LAE016 in the middle panel of 
FigureEl An emission line is clearly detected at A = 4973A, indicating a redshift of z = 3.091. The NB3640 detection 
is weak and offset by 0'.'8, but covered by the slit. As no additional emission features are detected in the spectrum, 
we retain LAE016 as a possible LyC-leaking galaxy 

LAE018 

The bottom panel of Figure [7] shows the HST / ACS image of LAE018 as well as the shallow and deep-mask spectra. 
An emission line at A = 4976A is detected in both spectra, indicating z = 3.093. The LyA flux is centered on a 
compact source to the north east, while the NB3640 flux extends over more diffuse non-ionizing UV flux to the south 
and west. Both slits cover the NB3640 flux. However, no other emission features were detected in either spectrum. 
Thus, we retain LAE018 as a possible LyC-leaking galaxy. 

LAE021 

As we do not have HST imaging of LAE021, we show our BV image and shallow- mask spectrum in the top panel 
of Figure [5J An emission line is clearly detected at A = 4948 A, indicating a redshift of z — 3.070. No other emission 
features are detected and we retain LAE021 as a possible LyC-leaking galaxy. We note, however, that the NB3640 
flux is offset by 1'.'4 from the LyA flux, and is not covered by the shallow-mask slit. As with all of our retained 
NB3640 detections, we account for the possibility the NB3640 flux is due to a foreground interloper in our Monte 
Carlo simulation. 

LAE028 

We obtained spectra of LAE028 on both shallow and deep masks. The middle panel of Figure [5] shows the H ST j ACS 
image. There is an offset of ~ 0'.'9 between the LyA and NB3640 flux centroids, each of which is centered on an 
individual clump. An emission line is detected at A ~ 4973A in both spectra indicating a redshift of z = 3.088. No 
additional emission or absorption features are detected in either spectra. We therefore retain LAE028 as a possible 
LyC-leaking galaxy. 

LAE038 

We obtained spectra of LAE038 in the deep mask as well as two of the shallow masks. The H ST / ACS image, the 
deep-mask spectrum and one of the shallow-mask spectra are shown in the bottom panel of Figure [5J An emission line 
at A ~ 4983A, indicating z — 3.099, is detected in each of the spectra. In all three of the spectra, the Lya emission line 
is offset from the continuum by ~ 0'.'5. The directions and magnitude of the offsets are consistent with the continuum 
emanating from the location of NB3640 flux and the Lya emission from the LyA flux. No additional lines are detected 
in any of the three spectra. We therefore retain LAE038 as a possible LyC-leaking galaxy. 

LAE039 

We show the BV image and shallow-mask spectrum of LAE039 in the top panel of Figure [9] We identify an emission 
line at A = 4978A, indicating z = 3.095. No additional lines are detected. We therefore retain LAE039 as a possible 
LyC-leaking galaxy. 

LAE041 

We obtained spectra of LAE041 in the deep mask as well as two of the shallow masks. The H ST / ACS image, the 
deep-mask spectrum and one of the shallow- mask spectra are shown in the middle panel of Figure [9l We identify an 
emission line at A ~ 4983A, indicating z = 3.066. No additional lines are detected in any of the three spectra. We 
therefore retain LAE041 as a possible LyC-leaking galaxy. 

LAE046 

We show our BV image and shallow-mask spectrum of LAE046 in the bottom panel of Figure [9l An emission line 
is detected at A = 4984A, indicating a redshift of z = 3.100. No other emission features are detected and we retain 
LAE046 as a possible LyC-leaking galaxy. We note, however, that the NB3640 flux is offset by 1'.'8 from the LyA 
flux, and is not covered by the shallow-mask slit. As with all of our retained NB3640 detections, we account for the 
possibility the the NB3640 flux is due to a foreground interloper in our Monte Carlo simulation. 
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LAE04S 

We show the H ST/ ACS image and shallow-mask spectrum of LAE048 in the top panel of Figure [TUJ An emission 
line is detected at A = 49 77 A, indicating a redshift of z — 3.094. No other emission features are detected and we 
retain LAE048 as a possible LyC-leaking galaxy. We note, however, that the NB3640 flux is offset by l'/7 from the 
LyA flux, and is not covered by the shallow-mask slit. As with all of our retained NB3640 detections, we account for 
the possibility the the NB3640 flux is due to a foreground interloper in our Monte Carlo simulation. 

LAE051 

We show the HST/ ACS imaging and shallow-mask spectrum of LAE051 in the middle panel of Figure [TU] An 
emission line is clearly detected at A = 4977A, indicating a redshift of z — 3.094. The shallow-mask slit also covers 
the position of the NB3640 detection. No additional emission features are detected in the spectrum. Thus, we retain 
LAE051 as a possible LyC-leaking galaxy 

LAE053 

We obtained spectra of LAE053 in one of the shallow masks and the deep mask. The bottom panel of Figure ITUl shows 
the H ST I ACS image and extracted spectra. There is an offset of ~ 0'.'9 between the LyA and NB3640 flux centroids, 
each of which is centered on an individual clump. Each clump and corresponding LyA or NB3640 flux is covered 
by only one of the masks. A relatively broad emission line is detected at A ~ 4973A in the shallow-mask spectrum, 
indicating a redshift of z = 3.088. An absorption feature is detected in the deep mask at A ~ 4967A indicating Lya 
in absorption at z = 3.085. This absorption feature suggests that the NB3640 source is also at z ~ 3.09. No other 
emission or absorption features are detected in either spectra. We therefore retain LAE053 as a possible LyC-leaking 
galaxy. 

LAE064 

We show the H ST/ ACS imaging and the combined shallow- mask spectrum of LAE064 in the top panel of Figure ITTl 
An emission line is clearly detected at A = 4994A, indicating a redshift of z = 3.108. Both slits also cover the position 
of the NB3640 detection, which is offset by l'.'O from the LyA flux. The extraction apertures for the spectra shown were 
centered on the Lya line, which was offset from faint continuum emission consistent in magnitude (~1") and direction 
from the position of the NB3640 flux. No additional emission features are detected in either (one- or two-dimensional) 
spectra. Thus, we retain LAE064 as a possible LyC-leaking galaxy. 

LAE069 

We show our BV image and shallow- mask spectrum of LAE069 in the middle panel of Figure [TT] An emission line 
is detected at A = 4953A, indicating a redshift of z — 3.074. No other emission features are detected and we retain 
LAE069 as a possible LyC-leaking galaxy. We note, however, that the NB3640 flux is offset from the LyA flux by 0'.'9, 
and is not covered by the shallow- mask slit. 

LAE074 

The bottom panel of Figure ITTl shows the HST / ACS image of LAE074 as well as the shallow and deep-mask spectra. 
An emission line at A = 4990A is detected in the shallow-mask spectrum, indicating z — 3.093. While the Lya emission 
line is not detected in the deep-mask spectrum, the slit covers only part of the diffuse LyA flux flux. No other emission 
features were detected in either spectrum. Thus, we retain LAE074 as a possible LyC-leaking galaxy. 

LAE081 

The top of Figure [12] shows the HST/ ACS image of LAE081 as well as the shallow and deep-mask spectra. An 
emission line at A = 4989A is detected in both spectra, indicating z = 3.104. Both slits cover the bulk of the LyA 
and NB3640 flux. No other emission features were detected in either spectrum. Thus, we retain LAE081 as a possible 
LyC-leaking galaxy. 

LAE083 

The bottom of Figure [T2"l shows the HST/ ACS image of LAE083 as well as the shallow-mask spectrum. An emission 
line at A = 4942A is detected, indicating z = 3.065. The slit also covers the bulk of the NB3640 flux. No other 
emission features were detected. Thus, we retain LAE083 as a possible LyC-leaking galaxy. 
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The z > 3.055 LBG sample. 
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D17 


22 


17 


18 


86 


00 


18 


17 





3 


090 


3 


070 


24 


29 


27 


00 


0' 


9 


12.2 ± 5.1 


M2 


22 


17 


33 


51 


00 


11 


10 


8 


3 


388 


3 


386 


25 


19 


25 


99 


0' 


8 


2.1 ± 0.4 


M5 


22 


17 


26 


07 


00 


11 


33 


1 


3 


327 






25 


71 


26. 


n e 


0' 


4 




M6 


22 


17 


28 


13 


00 


11 


40 


5 


3 


180 






25 


65 










> 4.6 


M8 


22 


17 


25 


10 


00 


11 


56 


8 


3 


064 


3 


062 


24 


56 










> 12.5 


M10 


22 


17 


26 


80 


00 


12 


21 


3 


3 


103 


3 


095 


24 


50 










> 13.1 


Mil 


22 


17 


31 


77 


00 


12 


51 


3 


3 


107 


3 


103 


25 


22 










> 6.8 


M13 


22 


17 


31 


46 


00 


12 


55 


2 


3 


107 






25 


38 










> 5.9 


M14 


22 


17 


39 


05 


00 


13 


30 


1 


3 


091 






25 


20 










> 6.9 


M19 


22 


17 


36 


90 


00 


15 


00 


9 






3 


082 


24 


96 










> 8.6 


M20 


22 


17 


34 


40 


00 


15 


02 


9 


3 


109 


3 


096 


25 


08 










> 7.7 


M25 


22 


17 


31 


49 


00 


16 


31 


2 


3 


098 


3 


091 


24 


79 










> 10.1 


M28 


22 


17 


31 


66 


00 


16 


58 





3 


094 


3 


088 


24 


75 










> 10.4 


M29 


22 


17 


37 


10 


00 


17 


08 


8 






3 


228 


24 


93 


26 


76 


1' 





5.4 ± 1.3 


M31 


22 


17 


36 


87 


00 


17 


12 


4 


3 


099 






25 


70 










> 4.4 


M34 


22 


17 


33 


80 


00 


17 


57 


2 






3 


084 


25 


41 










> 5.7 


MD14 


22 


17 


37 


91 


00 


13 


13 


9 






3 


097 


24 


49 










> 13.4 


MD23 


22 


17 


28 


01 


00 


14 


29 


6 


3 


085 


3 


075 


24 


34 










> 15.2 


MD32 


22 


17 


23 


70 


00 


16 


01 


6 


3 


102 






25 


14 


25. 


51° 


0' 


4 




MD46 


22 


17 


27 


28 


00 


18 


09 


7 


3 


091 


3 


080 


23 


49 


25 


22 


1' 





4.9 ±0.7 


aug96M16 f 


22 


17 


30 


86 


00 


13 


10 


8 










24 


47 


25 


23 


0' 


7 


2.0 ±0.4 



a Lyct emission redshift. 

k Interstellar absorption redshift. 

c Spatial offset between the centroids of i?-band and NB3640 emission. 

^ Observed ratio and uncertainty in non-ionizing UV and LyC flux densities, inferred from the 
NB3640— R color. This value has not been corrected for cither contamination by foreground sources or 
IGM absorption. 

c Spectrum contains evidence for the presence of a foreground object. Thus, the NB3640 flux is likely 
contaminated by non-ionizing UV flux from the interloper. 

We were unable to determine the redshift of the LBG aug96M16, which in INestor et all d201lD was 
erroneously associated with a nearby galaxy with z em — 3.298. We have removed it from our sample 
but include it in this table for completeness. 



NESTOR ET AL. 



TABLE 2 

LBGS WITH NEW REDSHIFTS 2 < 3.055. 



ID 


RA 


Dec 




2 (. b 


R 




(J2000) 


(J2000) 








C8 


22:17:27.64 


00:11:59.1 


2.459 




26.42 


C19 


22:17:38.36 


00:14:16.4 


2.993 




24.84 


D5 


22:17:34.65 


00:12:33.6 


2.699 


2.690 


25.24 


D10 


22:17:38.92 


00:14:32.1 




2.605 


24.61 


M3 


22:17:35.81 


00:11:16.5 


2.876 




25.27 


M16 


22:17:24.66 


00:14:07.8 


3.021 


3.018 


24.88 


M35 


22:17:32.67 


00:18:05.0 


2.938 




25.42 


MD11 


22:17:34.94 


00:13:22.7 




2.857 


24.79 


MD42 


22:17:35.83 


00:17:19.8 




2.792 


25.37 


MD45 


22:17:25.36 


00:18:04.2 


3.025 


3.022 


24.39 



a Lya emission redshift. 

b Interstellar absorption redshift. 



LYMAN-CONTINUUM EMISSIVITY AT z = 3.1 



TABLE 3 

Summary of LAE Data 



ID RA Dec R 

(J2000) (J2000) 



nm d.c 
UU1 ' 


09 
zz 


17.Qi) An 

1 ( .OZ.4U 


U 


1 1 


Q/l 1 
o4. 1 


oq no 


002 d '° 


22 


17:38.90 





11 


01.8 


24.28 


003 


22 


17:24.79 





17 


17.4 


24.42 


004 d -° 


22 


17:28.01 





14 


30.0 


24.34 


005 


22 


17:35.86 





15 


59.4 


25.65 


006 


22 


17:24.80 





11 


16.8 


> 27.0 


007 


22 


17:27.78 





17 


36.9 


> 27.0 


008 d - e 


22 


17:21.11 





15 


28.0 


24.87 


009 d 


22 


17:28.29 





12 


12.3 


25.84 


010 


22 


17:20.38 





18 


04.2 


25.77 


Oil 


22 


17:33.85 





12 


14.9 


26.07 


012 d ' e 


22 


17:31.69 





16 


57.6 


24.75 


013 


22 


17:27.18 





16 


21.7 


25.98 


014 


22 


17:19.25 





14 


50.9 


25.82 


015 


22 


17:21.84 





12 


12.7 


> 27.0 


016 


22 


17:35.61 





18 


00.2 


26.24 


017 


22 


17:25.40 





17 


16.8 


26.22 


018 


22 


17:39.01 





17 


26.4 


26.25 


019 


22 


17:26.15 





13 


20.1 


25.70 


020 


22 


17:37.33 





16 


31.4 


25.45 


021 


22 


17:18.77 





15 


18.1 


> 27.0 


022 


22 


17:19.68 





11 


49.4 


26.11 


023 


22 


17:31.73 





16 


06.9 


24.91 


024 


22 


17:34.17 





16 


09.7 


26.73 


025 


22 


17:36.74 





16 


28.8 


25.54 


026 


22 


17:18.96 





12 


00.8 


26.59 


027 


22 


17:24.94 





17 


17.3 


26.25 


028 


22 


17:31.80 





17 


17.9 


25.50 


029 d 


22 


17:31.49 





12 


55.0 


25.38 


030 


22 


17:21.75 





11 


38.8 


> 27.0 


031 


22 


17:33.63 





17 


15.1 


26.37 


032 


22 


17:26.61 





13 


18.1 


26.61 


033 


22 


17:37.50 





14 


08.3 


26.62 


034 


22 


17:23.41 





16 


35.4 


25.42 


035 


22 


17:27.03 





13 


13.2 


> 27.0 


036 


22 


17:22.25 





11 


55.1 


> 27.0 


037 


22 


17:20.96 





18 


07.3 


25.61 


038 


22 


17:34.77 





15 


41.3 


26.17 


039 


22 


17:24.08 





11 


31.7 


26.48 


040 


22 


17:31.93 





13 


08.5 


> 27.0 


041 


22 


17:24.54 





15 


06.7 


25.97 


042 


22 


17:21.50 





17 


04.7 


25.50 


043 


22 


17:21.65 





12 


23.4 


26.24 


044 


22 


17:36.41 





12 


51.0 


> 27.0 


045 


22 


17:35.97 





16 


30.2 


> 27.0 


046 


22 


17:21.47 





14 


54.6 


>27.0 


047 


22 


17:36.05 





15 


06.9 


> 27.0 


048 


22 


17:27.37 





16 


51.5 


26.73 


049 


22 


17:39.29 





16 


10.5 


> 27.0 


050 


22 


17:24.56 





15 


56.8 


> 27.0 


051 


22 


17:33.72 





15 


04.9 


26.26 


052 


22 


17:36.84 





13 


17.2 


26.63 


053 


22 


17:34.70 





16 


33.4 


26.53 


054 


22 


17:39.05 





11 


33.9 


26.25 


055 


22 


17:35.80 





11 


50.0 


26.67 



NB3640 






Fuv b 

F LyC obs 




> 27.3 






> 22.4 


3.075 


> 27.3 






> 16.2 


3.076 


24.74 f 


0'.'3 


/ . / 6 




3.097 


07 Q 

> Z / .0 






> 10. z 


q nn9 


— 97 Q 
> Z / .0 






> 4.0 


o.uyo 


\ 97 *3 
> Z 1 .0 








q n7i 
o.U / 1 


^» 97 1 
.> z / .0 








o.uoy 


\ 97 3 
> Z I .0 






\ Q Q 


o.U/D 


> 27.3 






> 3.8 


3.071 


26.74 


n"o 

. 3 


n't o 

. 3 


2.4 ±1.1 


3.096 


> 27.3 






> 3.1 


3.102 


> 27.3 






> 10.5 


3.094 


> 27.3 






> 3.4 


3.095 


> 27.3 






> 3.9 


3.067 


> 27.3 








3.067 


26.91 


nff n 
0.9 


n// o 

0.8 


1.8 ± 0.9 


3.091 


\ 97 Q 

> At .0 






\ O 7 
> z. i 


o. 1UO 


zo.oy 


n"i 


i "n 
1 . u 


U.D ± U.Z 


o.uyo 


26.24 f 


Cft A 

0.4 


nil ( i 

0.9 




3.101 


> 27.3 








• ■ • S 


27.17 


A ft O 

1.3 


A II A 

1.4 


< 1.8 


3.070 


> 27.3 






> 3.0 


3.066 


— 97 Q 
> Z / .0 






> y.i 


o. 1U1 


> 27.3 






\ 1 7 

> 1.7 


3.096 


25.85 f 


0'.'3 


1"2 




3.091 


> 27.3 






> 1.9 


3.092 


> 27.3 






> 2.6 


3.101 


26.71 


0.3 


nil n 
0.9 


3.1 ± 1.3 


3.088 


> 27.3 






> 5.9 


3.107 


> 27.3 








3.096 


> 27.3 






> 2.3 


3.092 


> 27.3 






> 1.9 


3.101 


> 27.3 






> 1.9 


3.076 


25.76 h 


0'.'5 


o'.'o 


1.4 ±0.4 


3.044 


> 27.3 








3.093 


> 27.3 








3.084 


> 27.3 






> 4.7 


3.089 


25.82 


O'.'l 


0'.'7 


0.7± 0.3 


3.099 


26.77 


0'.'4 


0'.'8 


1.3 ±0.7 


3.095 


> 27.3 








3.105 


25.94 


0'.'4 


0'.'5 


1.0 ±0.4 


3.066 


> 27.3 






> 5.2 


3.072 


> 27.3 






> 2.7 


3.071 


> 27.3 








3.059 


> 27.3 








3.094 


26.43 


l'/8 


l'/8 


< 0.8 


3.100 


> 27.3 








3.096 


26.00 


2'/2 


l'/7 


0.5 ±0.2 


3.094 


> 27.3 








3.094 


> 27.3 








3.082 


27.21 


0'/3 


0'/7 


2.4 ± 1.3 


3.094 


> 27.3 






> 1.8 


3.102 


26.98 


0"6 


0"9 


1.5 ±0.8 


3.090 


> 27.3 






> 2.6 


3.096 


> 27.3 






> 1.8 


3.072 



NESTOR ET AL. 



TABLE 3 

Summary of LAE Data ( Continued) 



ID 


(J2000) 


JJCC 

(J2000) 


p 


IN IJOO'-tU 


A a 




^(7 V b 


~ c 


056 


22:17:22.42 


0:17:20.7 


26.86 


> 27.3 








> 1.5 


3.073 


057 


22:17:2&A0 


0:10:58.3 


26.84 


> 27.3 








> 1.5 


3.070 


058 


22:17:19.61 


0:15:38.4 


> 27.0 


> 27.3 










... 1 


059 


22:17:24.98 


0:12:30.0 


25.31 


> 27.3 








> 6.2 


3.096 


060 


22:17:28.19 


0:11:17.1 


26.61 


> 27.3 








> 1.9 


3.065 


061 


22:17:34.10 


0:15:40.2 


> 27.0 


> 27.3 










3.101 


062 


22:17:22.87 


0:14:41.7 


26.53 


> 27.3 








> 2.0 


3.057 


063 


22:17:23.32 


0:15:52.9 


26.55 


> 27.3 








> 2.0 


3.099 


064 


22:17:35.42 


0:12:14.6 


26.05 


26.74 


r\// 1 
. 1 


1' 


'0 


1.9 ± 0.9 


3.108 


065 


22:17:28.15 


0:14:oo.4 


26.91 


> 27.3 








> 1.4 


3.102 


066 


00.1 *7. on Oi^ 
22: 1 / :2U.oo 


0:15:11.8 


26.64 


> 27.3 








> 1.8 


3.066 


067 


00.1 T.oa on 
22: 1 / :oo.2o 


0:13:11.7 


26.40 


> 27.3 








> 2.3 


3.106 


068 


OO 17 1 O OT 

22:17:18. o7 


0:17:26.1 


> 27.0 


> 27.3 








■ • • 


3.072 


069 


00.1 7.1 o o^; 
22:1 /:18.yo 


0:11:12.0 


24.62 


27.22 


U . 2 


0' 


'9 


10.9 ± 5.1 


3.074 


u / u 


OO.I 7.QO OQ 

zz. i / .oy.zo 


U. 14.UU.Z 


9^ Ofi 

Zo.yo 


> Z / .O 








> 0.4 


Q AOA 

o.uyu 


071 
U i 1 


99-17-91 R1 


fi-1 9-90 ^ 

u. ±z. zu.o 


9fi 77 
ZD. 1 1 


^ 97 1 
.> z / .O 








S> l.D 


1 f!79 

o.u / z 


072 


99-1 7-^1 94 


0:17:32.1 


27.00 


> 27.3 








^ 1 ^ 


3.084 


073 


22:17:39.12 


0:17:11.7 


26.35 


> 27.3 








> 2.4 


3.083 


074 


22:17:36.47 


0:12:54.8 


26.16 


25.52 


O'/l 


u 


'7 
1 


0.6 ± 0.2 


3.105 


075 


22:17:22.97 


0:11:25.8 


> 27.0 


> 27.3 












076 


22:17:20.67 


0:15:13.2 


> 27.0 


> 27.3 










3.066 


077 


22:17:37.95 


0:11:01.3 


26.01 


26.36 


0"2 


1 


1 






078 


22:17:37.68 


0:16:48.3 


25.95 


> 27.3 








> 3.5 


3.090 


079 


22:17:34.68 


0:11:10.5 


> 27.0 


> 27.3 










3.103 


080 


22:17:35.95 


0:13:43.3 


> 27.0 


> 27.3 










3.104 


081 


22:17:29.22 


0:14:48.7 


>27.0 


26.79 


O'.'6-i 


u 


' 7 


< 1.2 


3.104 


082 


22:17:35.44 


0:16:47.6 


> 27.0 


> 27.3 










3.087 


083 


22:17:28.46 


0:12:08.9 


26.46 


26.84 


0"4 


u 


D 


1.4 ±0.7 


3.065 


084 


22:17:19.90 


0:15:14.9 


>27.0 


26.50 


O'.'l 


0' 


'8 




..." 


085 


22:17:30.86 


0:14:38.2 


26.94 


> 27.3 










... 1 


086 


22:17:28.42 


0:13:42.8 


> 27.0 


> 27.3 










3.100 


087 


22:17:37.07 


0:13:21.5 


>27.0 


27.26 


O'.'IJ 


1' 


'6 




. . . s 


088 


22:17:38.45 


0-13-18 3 


> 27.0 


> 27.3 










i 


089 


22-1 7-38 54 


0-1 5-22 5 


26.59 


> 27.3 










i 


090 


22:17:18.25 


0:14:06.4 


> 27.0 


> 27.3 










3.095 


091 


22:17:36.14 


0:15:40.7 


> 27.0 


> 27.3 










3.095 




99-1 7-9^ Q7 
zz. i ( .zo.y < 


n-1 ^-97 & 


^> Z 1 .u 


\ 97 q 
^> z / .o 










i 


uyo 


99-1 7-97 AS 

ZZ. l/.Zl ."4:0 


fl-1 ^- ^7 ^ 
U. lO. I . 


9fi IX. 
ZO.OO 


- , 97 o 
s> Z i .o 








^ Z.4 


o. 1UD 


uy4 


99-1 7.QQ 1 /I 

zz. i / .oy. 14 


n. i 7-nn ^ 


97 n 

> Z / .U 


\ 97 "} 
> Z / .o 












uyo 


99-1 l-'XI 1 Q 
zz. i/.oi. if 


n-1 3-9& n 


9fi 1 7 
ZU. ± ( 


- , 97 q 








^> Z.O 


o.uy4 


uyo 


99-1 7-^S 

zz. i ( .oo.yo 


n-1 1 -^7 a 


9fi 

ZD.D4 


9^ 


n"fi 

U . D 


0' 


'6 




g 


uy / 


99-1 7-97 1 1 
ZZ. li.Zl.ll 


U.14.U0. i 


> Z / .u 


\ 97 Q 










O.UD4 


uyo 


99-1 7-9/1 fl"l 
ZZ. 1 ( .Z4.U1 


fl.1Q.1Q c 


\ 97 n 
^ Z / .u 


\ 97 Q 
> Z ( .O 










g 


uyy 


99-1 7-^fi Afi 
ZZ. 1 ( .OU. i ±U 


n-1 ^ nn ^ 

U. -LO.UU.O 


, 97 n 
z f .u 


^ 97 ^ 
^ z / .o 










*\ 1 09 
0. 1UZ 


100 


22:17:30.61 


0:18:11.6 


> 27.0 


> 27.3 












101 


22:17:25.33 


0:17:22.5 


26.89 


26.58 


0"2 


1' 


'4 




. . . s 


102 


22:17:24.00 


0:16:27.6 


25.91 


26.05 


O'.'l 


0' 


'2 




. . . s 


103 


22:17:19.40 


0:15:26.1 


> 27.0 


> 27.3 










3.099 


104 


22:17:37.66 


0:12:55.5 


> 27.0 


> 27.3 










3.097 


105 


22:17:35.46 


0:12:23.9 


26.81 


> 27.3 








> 1.6 


3.104 


106 


22:17:22.86 


0:17:57.8 


26.16 


> 27.3 












107 


22:17:20.96 


0:14:46.7 


26.57 


> 27.3 








> 2.0 


3.075 


108 


22:17:24.78 


0:17:40.4 


26.84 


> 27.3 








> 1.5 


3.100 


109 


22:17:23.98 


0:17:57.8 


26.77 


> 27.3 










. . . s 


110 


22:17:19.50 


0:15:57.6 


> 27.0 


> 27.3 










3.074 



LYMAN-CONTINUUM EMISSIVITY AT z = 3.1 



TABLE 3 

Summary of LAE Data ( Continued) 



ID 


RA 
(J2000) 


Dec 
(J2000) 


R 


NB3640 




A LyA 


Fjjy b 
F LyC l s 




111 


22 


17 


31 


14 





16 


42 


9 


> 27.0 


> 


27.3 








3 


096 


112 


22 


17 


32 


72 





15 


54 


2 


> 27.0 


> 


27.3 








3 


096 


113 


22 


17 


24 


80 





13 


26 


9 


> 27.0 


> 


27.3 








3 


069 


114 


22 


17 


34 


50 





14 


20 





> 27.0 


> 


27.3 








3 


103 


11 r 

llo 


22 


17 


33 


46 





17 


01 


2 


> z/.U 


> 


27.3 








3 


093 


110 


22 


17 


28 


00 


o 


12 


14 


2 


07 n 


> 


27.3 








3 


096 


117 
11 I 


22 


17 


39 


08 





12 


01 


9 


\ 07 n 
> Zi .U 


> 


27.3 




"' 




3 


098 


118 


22 


17 


35 


28 





10 


59 


9 


26.20 


26.11 


. 4 


1 . 4 






. . g 


119 


22 


17 


25 


63 





12 


47 


8 


> 27.0 


> 


27.3 








3 


071 


120 


22 


17 


26 


76 





10 


59 


8 


> 27.0 


> 


27.3 








3 


079 


121 


22 


17 


26 


44 





15 


27 


5 


> 27.0 


> 


27.3 








3 


099 


122 


22 


17 


38 


19 





14 


03 


7 


> 27.0 


> 


27.3 








3 


106 


123 


22 


17 


35 


06 





17 


26 





> 27.0 


> 


27.3 












124 


22 


17 


22 


80 





17 


48 


7 


> 27.0 


> 


27.3 










. . s 


125 


22 


17 


38 


02 





14 


03 


6 


26.79 


> 


27.3 






> 1.6 


3 


101 


126 


22 


17 


19 


53 





16 


48 


2 


> 27.0 


> 


27.3 












127 


22 


17 


36 


91 





11 


27 


1 


> 27.0 


> 


27.3 








3 


107 


128 


22 


17 


23 


43 





16 


07 


4 


> 27.0 


> 


27.3 








3 


095 


129 


22 


17 


22 


28 





10 


57 


9 


> 27.0 


> 


27.3 








3 


111 


130 


22 


17 


32 


84 





16 


48 


8 


> 27.0 




27.3 








3 


092 



a Spatial offset between the centroids of R or Lya and NB3640 emission. 

b Observed ratio and uncertainty in non-ionizing UV and LyC flux densities, inferred from the NB3640— R color. 
This value has not been corrected for either contamination by foreground sources or IGM absorption. 
c Spectroscopic redshift based on the observed wavelength of Lya emission. 

d Objects that arc also identified as LBGs. 001 is LBG D3; 002 is LBG C4; 004 is LBG MD23; 008 is LBG C28; 

009 is LBG C9; 012 is LBG M28; 029 is LBG M13. 

c Objects with previously-known spectroscopic redshifts. 

^ Spectrum contains evidence for the presence of a foreground galaxy. Thus, we do not consider the NB3640 
detection to be LyC flux. 

g Objects with LRIS spectroscopic observations but without confirmed redshifts. 

k The redshift of this LAE is such that some non-ionizing UV flux will contribute to the NB3640 detection, and 

therefore we exclude it from our LyC detection sample. 

1 Objects without LRIS spectroscopic observations. 

J Offset determined from the ccntroid of the BV detection. 
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TABLE 4 

Sample average colors and UV to LyC flux density ratios. 





LBGs 




LAEs 




Correction 


(NB3640) - (i?) a 




(NB3640) - (R) a 






Full Ensembles 






none 

contamination 11 

IGM + contamination 


3 74 +0.66 
°- ' -0.42 

4 45+1.16 


oi r+26.4 
OJ "°-1Q.O 

eo.s™ 

18.0t 7 4 4 8 


9 -, ,+0.36 
^' ±o -0.29 

2 66+ 0,56 
z - DD -0.39 

1 42+ - 57 
^ -0.40 


7 1+2.8 
'• x -1.7 

n.6±I;l 

o y+2.5 


Full Ensembles, 


maximum allowed flux for NB3640 non-detections' 




none 

contamination" 1 

IGM + contamination 


3 40+ 43 

3 84+ - 48 
9 rq+0.51 
Z - JO -0.37 


23.0±ft° 
34.4+^3 

10.2±f;J 


1 90+ 029 

i.»u_ 25 
9 9 o+0.37 
^S_0.29 
1 04+ 038 


5.7+ 1 ' 7 

2 6+ 1 1 
Z - D -0.6 




NB3640 detections only 






none 

contamination d 

IGM + contamination 


i 77 +0.67 
'-0.63 

1 8Q+ - 75 
1 -° a -0.66 
57+O.77 
U -°'-0.68 


5 1+4.4 
°- 1 -2.3 
5 7+5.7 
°-'-2.6 

1 7+ 1 - 7 


29+ - 41 
u,ZB -0.39 
n OO+0.46 
u -°°-0.42 
-0 91+ - 47 


i q+0.6 
1 - J -0.4 

^-o^ 

4+ ' 2 
u -*-o.i 



a Color determined from average NB3640 and i?-band fluxes. Uncertainties include individual 
flux and sample uncertainties. 

k Ratio and uncertainty in non-ionizing UV and LyC flux density inferred from (NB3640) — (R) 
color. 



c Color and flux-density ratio determined assuming NB3640 non-detections contribute zero 
LyC flux. 

^ Color and flux-density ratio after statistically correcting sample for foreground contamination 
of NB3640 fluxes. 

e Color and flux-density ratio after correcting sample for foreground contamination and IGM 
absorption of NB3640 fluxes. 

^ Color and flux-density ratio determined assuming NB3640 non-detections contribute the max- 
imum possible LyC flux consistent with our limits from stacking analysis. 



TABLE 5 

Contributions to the ionizing backgrounds. 



LF a rj h Magnitude range c € LyC d 



(i) 


LBG 


18.0± 7 4 4 8 


M A b < -20.0 


c n+3.6 
°' z -3.4 


(ii) 


LAE 


q 7+2.5 
°-'-l.l 


-20 < M AB < -18.3 


fi 7+ 2 - 8 


(iii) 


LBG 


n 7+2.5 
^-'-l.l 


-20 < M A b < -18.3 


27 0+ 11 ' 4 
z '- u -10.9 


(iv) 


LBG 


18.0l 7 4 4 8 


M AB < -18.3 


10.7+ 7 f 


(v) 


LAE 


q 7+2.5 

"•'-1.1 


M A b < -18.3 


8-6+^ 




Total (lum.-dep.) e 




M A b < -18.3 


o 9 9 +12.0 
■"-11.4 




Total (LAE-dep.) f 




M A b < -18.3 


I6.8j£j 



a Luminosity function parameters used in Equations [3V-f5l LBG from Rcddv ct al. (2008); 
or LAE from lOuchi et all J2008fl scaled to include LAEs having REW> 20A. 



k Sample average flux-density ratio used in Equations [3T-f5l 
c Magnitude range over which the first moment of the luminosity function is determined. 
Mab — —20.0 and —18.3 correspond to QA6L^ BG and 0.1L^ SG , respectively. 
^ Comoving specific emissivity of ionizing radiation in units of 10 24 ergs s — 1 Hz -1 Mpc -3 . 
c Totals for our luminosity-dependent r} model, determined by summing rows (i) and (iii). 

Total for our LAE-dcpcndcnt r/ model, determined by summing 0.77x row (iv) and row 
(v). 



LYMAN-CONTINUUM EMISSIVITY AT z = 3.1 



TABLE 6 

Estimates of the intergalactic hydrogen 
photoionization rate at z — 3.1. 



Reference 



r (xio 



-12 c 



- 1 ) 



Meiksin fe White (20041 
Bolton fc Haehnelt (2007) 
Faucher-Giguere et al. (2008 1 
This work, luminosity-dependent r) 
This work, LAE-depcndent rj 



a +0.l4 
-0.12 
q+0.34 
-0.26 

0.63 ±0.08 

7+1-0 
-0.9 
1+0.6 
-0.6 



0.88Z 



2.72 
1.4-* 



TABLE 7 

Model dependent inferred values for /, 



LyC 









BC03 a 






BPASS b 




Age (yr) c 


z d 


y-stars c 


f LyC,LBG f 
J esc 


J esc 


YjStars c 


rLyC,LBG f 
Jesc 


f -LyC,LAE g 
J esc 








fUV 
J esc 


= 0.2 








10 (> 


0.004 


1.98 


0.02 


0.11 


1.33 


0.01 


0.07 




0.020 


1.90 


0.02 


0.10 


1.45 


0.02 


0.08 


10 7 


0.004 


3.59 


0.04 


0.19 


2.10 


0.02 


0.11 




0.020 


4.20 


0.05 


0.23 


2.70 


0.03 


0.15 


10 s 


0.004 


6.17 


0.07 


0.33 


3.16 


0.04 


0.17 




0.020 


6.38 


0.07 


0.34 


4.43 


0.05 


0.24 








fUV 
J esc 


= 0.3 








10 (> 


0.004 


1.98 


0.03 


0.16 


1.33 


0.02 


0.11 




0.020 


1.90 


0.03 


0.15 


1.45 


0.02 


0.12 


10 7 


0.004 


3.59 


0.06 


0.29 


2.10 


0.04 


0.17 




0.020 


4.20 


0.07 


0.34 


2.70 


0.05 


0.22 


10* 


0.004 


6.17 


0.10 


0.50 


3.16 


0.05 


0.26 




0.020 


6.38 


0.11 


0.52 


4.43 


0.07 


0.36 








r uv 

J esc 


= 1.0 








10 8 


0.004 


1.98 


0.11 


0.54 


1.33 


0.07 


0.36 




0.020 


1.90 


0.11 


0.97 


1.45 


0.08 


0.39 


10 7 


0.004 


3.59 


0.20 


> 1 


2.10 


0.12 


0.57 




0.020 


4.20 


0.23 


> 1 


2.70 


0.15 


0.73 


10 8 


0.004 


6.17 


0.34 


> 1 


3.16 


0.18 


0.85 




0.020 


6.38 


0.35 


> 1 


4.43 


0.25 


> 1 



Stanwav (2009). 



k Stellar population synthesis models of Eldridec . 
c Time since turn on of constant star formation. 
^ Metallicity of stellar population. 

c Ratio of intrinsic i^ieoo to -Fnb3640 predicted by model. 

^ Resulting LyC escape fraction for the given UV escape fraction, using tilbg- 
s Resulting LyC escape fraction for the given UV escape fraction, using tilae- 



